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TECHNICAL  NOTE  2lt.09 

SUMMARY  OF  METHODS  FOR  CALCULATING  DYNAMIC  LATERAL 
STABILITY  AND  RESPONSE  AND  FOR  ESTIMATING 
LATERAL  STABILITY  DERIVATIVES 
By  John  P.  Campbell  and.  Marion  0.  McKinney 

SUMMARY 


A  sunmary  of  methods  for  making  dynamic  lateral  stability  and 
response  calculations  and  for  estimating  the  aerodynamic  stability 
derivatives  required  for  use  in  these  calculations  is  presented.  The 
processes  of  performing  calculations  of  the  time  histories  of  lateral 
motions,  of  the  period  and  damping  of  these  motions,  and  of  the  lateral 
stability  boundaries  are  presented  as  a  series  of  simple  straightforward 
steps.  Existing  methods  for  estimating  the  stability  derivatives  are 
summarized  and,  in  some  cases,  simple  new  empirical  formulas  are  pre¬ 
sented.  Reference  is  also  made  to  reports  presenting  experimental  data 
that  should  be  useful  in  making  estimates  of  the  derivatives.  Detailed 
estimation  methods  are  presented  for  low-subsonic-speed  conditions  but 
only  a  brief  discussion  and  a  list  of  references  are  given  for  transonlc- 
and  supersonic-speed  conditions. 


INTRODUCTION 


Dynamic  lateral  stability  has  not  received  widespread  attention  in 
the  past  because  it  has  not  generally  been  a  serious  problem  in  the 
design  of  airplanes.  Consideration  of  dynamic  lateral  stability  has 
recently  become  more  important,  however,  because  current  design  trends 
toward  the  use  of  low  aspect  ratio,  sweepback,  and  higher  wing  loading 
have,  in  many  cases,  led  to  unsatisfactory  dynamic  lateral  stability. 
Airplane  designers  are  therefore  finding  it  necessary  to  make  such  calcu¬ 
lations  in  connection  with  the  design  and  modification  of  airplanes.  In 
many  cases  these  calculations  are  difficult  to  perform  for  designers  who 
have  had  no  previous  experience  in  theoretical  stability  work  because 
most  of  the  published  theoretical  analyses  are  not  presented  in  a  forai 
that  is  especially  suited  to  the  computation  of  dynamic  stability.  The 
estimation  of  the  stability  derivatives  required  in  dynamic  stability 
calculations  has  also  been  found  to  be  difficult  in  many  cases.  Although 
theoretical  and  experimental  data  on  these  derivatives  have  appeared  in 
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numerous  publications^  no  single  publication  has  presented  methods  for 
estimating  the  derivatives  for  all  types  of  airplanes. 

One  approach  to  a  presentation  of  methods  of  calculating  stability 
and  estimating  stability  derivatives  in  a  form  suitable  for  use  by 
designers  was  made  by  Zimmerman  in  reference  1.  Although  this  report 
has  proved  to  be  of  valuable  assistance  to  designers  in  maJ5:ing  dynamic 
stability  calculations,  recent  trends  in  airplane  design  have  caused  its 
usefulness  to  be  seriously  limited.  For  example,  the  equations  of  refer¬ 
ence  1  do  not  include  the  product-of-inertia  terms  which  have  been  shown 
by  recent  studies  to  be  very  important  in  some  cases.  (See  references  2 
and  3.)  Moreover,  the  calculation  of  the  time  histories  of  lateral 
motions,  one  type  of  calculation  that  has  been  the  subject  of  increasing 
interest  in  the  last  few  years  (references  4  to  7)^  is  not  covered  in 
reference  1.  The  methods  of  estimating  stability  derivatives  presented 
in  reference  1  are  also  limited  because  they  apply  only  to  airplanes 
having  unswept  wings  with  an  aspect  ratio  of  6  operating  at  speeds  at 
which  compressibility  effects  are  negligible.  The  purpose  of  the  present 
paper  is  to  extend  the  methods  of  reference  1  to  include  the  methods  of 
computation  which  are  of  current  interest  to  designers  and  to  include 
methods  of  estimating  derivatives  for  configurations  and  flight  conditions 
which  are  now  being  considered. 

This  paper  summarizes  and  reduces  to  simple  strai^t  for  ward  steps 
methods  for  computing  the  time  histories  of  lateral  motions,  the  period 
and  damping  of  these  motions,  and  the  lateral  stability  boundaries. 
Existing  methods  of  estimating  stability  derivatives  for  a  variety  of 
airplane  configurations  are  summarized  and,  in  some  cases,  simple  new 
empirical  formulas  are  presented.  Reference  is  also  made  to  reports 
presenting  experimental  data  that  should  be  useful  in  making  estimates 
of  these  derivatives. 


SYMBOLS 

All  forces  and  moments  are  referred  to  the  stability  system  of  axes 
which  is  defined  in  figure  1.  The  following  definitions  apply  to  the 
symbols  except  where  they  are  otherwise  defined: 

m  mass  of  airplane,  slugs 

S  wing  area,  square  feet 

c  wing  mean  chord,  feet  (b/A) 

b  wing  span,  feet 

y^  span  of  that  part  of  wing  that  has  tip  dihedral,  feet 
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tail  length  (distance  from  center  of  pressure  of  vertical 
tail  to  center  of  gravity,  measured  parallel  to  longi¬ 
tudinal  stability  axis;  values  of  I  must  be  calculated 
for  each  angle  of  attack),  feet 

average  fuselage  height  at  ving  root,  feet 

average  fuselage  width  at  wing  root,  feet 

vertical  distance  of  quarter  chord  of  wing  root  chord  from 
fuselage  center  line,  positive  downward,  feet 

nondimensional  time  parameter  based  on  span  (Vt/b) 

longitudinal  distance  rearward  from  airplane  center  of 
gravity  to  wing  aerodynamic  center,  feet 

longitudinal  distance  from  leading  edge  of  vertical  tail 
chord  to  horizontal  tail  aerodynamic  center,  feet 
(see  fig.  6) 

vertical  distance  from  horizontal  tail  to  base  of  vertical- 
tail,  feet  (see  fig.  6) 

height  of  center  of  pressure  of  vertical  tail  above  longi¬ 
tudinal  stability  axis;  values  of  z  must  be  calculated 
for  each  angle  of  attack,  feet 

aspect  ratio 

sweepback  of  wing  quarter-chord  line,  degrees 

taper  ratio  (Tip  chord/Root  chord);  also,  differential 
operator  in  Laplace  transform 

dihedral  angle,  degrees  (see  sketch  of  fig.  9) 

dihedral  angle  of  -  wing  tip,  degrees 

time,  seconds 

airspeed,-  feet  per  second 

radius  of  gyration  about  principal  longitudinal  axis  of 
inertia,  feet 

radius  of  gyration  about  principal  normal  axis  of  inertia, 
feet 
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radius  of  gyration  about  X  axis^  feet 


k^/b 

product- of- inertia  factor 


V 

angle  of  attack  of  principal  longitudinal  axis  of  inertia, 
degrees  (see  fig.  2) 

angle  of  climb,  degrees  (see  fig.  2) 

angle  of  attack  of  longitudinal  body  axis,  degrees 
(see  fig,  2) 

angle  between  principal  longitudinal  axis  of  inertia  and 
longitudinal  body  axis,  degrees  (see  fig,  2) 

air  density,  slugs  per. cubic  foot 

angle  of  bank,  radians 
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angle  of  yaw,  radians 
P  angle  of  sideslip,  radians 

p  rolling  velocity,  radians  per  second  (d0/dt) 

r  yawing  velocity,  radians  per  second  (dilf/dt) 

0o  initial  angle  of  bank,  radians 

to  initial  angle  of  yaw,  radians 

Pq  initial  angle  of  sideslip,  radians 

(D0)q  nondimensional  initial  rolling  velocity  (d0/d0) 

(^t)o  nondimensional  initial  yawing  velocity  (dt/dc) 

R  Routh's  discriminant  or  real  part  of  complex  root  R  +  li 

I  imaginary  part  of  complex  root  R  +  li 

A,B,C,D,E  coefficients  of  the  characteristic  biquadratic  equation 
•••  factors  of  the  B,  C,  and  D  coefficients 

roots  of  characteristic  biquadratic  equation 
D  differential  operator  (d/dc) 

P  period  of  the  lateral  oscillation,  seconds 

^1/2  "time  to  damp  to  one -half  amplitude,  seconds 

T  time  conversion  factor  (m/pSV) 

o  nondimensional  time  factor  (t/r) 

H  relative  density  factor  (m/pSb) 

^c  impressed  rolling  moment,  foot-pounds 

^c  Impressed  yawing  moment,  foot-pounds 

^c  impressed  lateral  force,  pounds 
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impressed  rolling-moment  coefficient 
impressed  yawing -moment  coefficient 
impressed  lateral -force  coefficient 
lift  coefficient  (Lift/qS) 
drag  coefficient  (Drag/qS) 

rolling -moment  coefficient  (Rolling  moment/qSb) 
yawing -moment  coefficient  (Yawing  moment/qSb) 
lateral-force  coefficient  (Lateral  force/qS) 
dynamic  pressure,  pounds  per  square  foot 
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V  =  iir 


incremen't  in  produced  by  lift  and  induced-drag 

forces 

increment  in  Crip  produced  by  drag  not  associated  with 
lift  ^ 


Subscripts 


horizontal  tail 


section  lift  curve  slope 


wing 


wing 


design 


fuselage 

used  to  designate  vertical  tail 

used  to  designate  design  under  consideration 

used  to  designate  design  for  which  force-test  data  are 
available 

experimental 


V-tail 


V-tail 


effective 


horizontal  tail 


NACA  TN  2409 


9 


CALCULATION  OF  LATERAL  STABILITY 

AND  RESPONSE 


Various  types  of  calculations  may  be  performed  to  indicate  in  some 
way  the  stability  of  an  airplane  or  the  response  to  gust  disturbances 
and  control  manipulations.  The  calculations  most  commonly  made  are  cal¬ 
culations  of  time  histories  of  disturbed  motions^  period  and  damping  of 
the  free  motions^  and  spiral  and  oscillatory  stability  boundaries  (lines 
of  neutral  damping  of  the  spiral  mode  and  of  the  lateral  oscillations). 
Step-by-step  procedures  for  performing  these  types  of  calculations  are 
explained  in  the  text  and  derivations  and  additional  pertinent  material 
are  presented  in  appendixes  A  to  D. 

The  period  and  damping  calculations  are  the  easiest  of  the  three 
types  to  perform.  For  this  reason^  and  because  the  dynamic  lateral 
stability  of  airplanes  is  at  present  specified  in  the  flying-qualities 
requirements  in  terms  of  the  period  and  damping  of  the  lateral  oscilla¬ 
tion^  period  and  damping  calculations  are  probably  the  most  commonly 
performed. 

Recent  dynamic  stability  work  has  indicated^  however^  that  the 
period  and  damping  characteristics  of  the  free  motions  of  an  airplane  . 
are  not  always,  a  sufficient  indication  of  whether  the  dynamic  behavior 
of  an  airplane  following  various  types  of  disturbances  will  be  con¬ 
sidered  satisfactory.  For  this  reason  the  calculation  of  time  histories 
of  the  motions  of  airplanes  is  becoming  more  common  despite  the  fact 
that  these  calculations  are  fairly  laborious.  The  increasing  use  of 
automatic  computing  machines  has  also  made  the  calculation  of  motions 
more  popular. 

For  many  years calculations  of  stability  boundaries  were  the  type 
of  calculation  most  commonly  performed.  In  recent  years ^  however,  sta¬ 
bility  boundaries  have  not  been  considered  to  give  an  adequate  indica¬ 
tion  of  stability.  Since  boundaries  are  useful  in  some  cases,  however, 
(for  example,  for  quick  approximation  of  the  effects  of  changes  in 
dihedral  and  tail  area)  the  methods  of  calculating  the  spiral  and  oscil¬ 
latory  stability  boundaries  are  described  herein.  Lines  of  constant 
period  and  damping  of  the  lateral  oscillation  are  related  to  stability 
boundaries  (lines  of  neutral  stability).  In  some  cases  these  lines  of 
constant  period  and  damping  may  prove  more  useful  than  boundaries.  Since 
no  extensive  use  has  been  made  of  lines  of  constant  period  and  damping, 
however,  the  methods  of  calculating  these  lines  (presented  in  refer¬ 
ences  8  and  9)  are  not  given  in  the  present  paper. 

The  equations  and  methods  of  calculation  presented  in  the  present 
paper  deal  specifically  with  the  inherent  motions  of  airplanes  for  the 
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case  of  three  degrees  of  freedom  (roll,  y&w,  and  sideslip)  and  linear 
stability  derivatives.  In  order  to  perform  similar  calculations  for 
cases  involving  additional  degrees  of  freedom,  nonlinear  derivatives, 
or  autopilots  with  time  lag,  special  equations  are  required.  The 
methods  and  equations  for  treating  these  cases  are  presented  in  refer¬ 
ences  10  to  18.  Additional  degrees  of  freedom  for  the  case  of  free 
controls  are  treated  in  references  l6  to  l8  and  for  the  case  of  fue 
sloshing  are  treated  in  reference  10.  The  use  of  nonlinear  derivatives 
in  stability  calculations  is  covered  in  reference  11.  Methods  of 
treating  the  effect  of  autopilots,  including  the  effect  of  time  lag  in 
the  autopilot  are  presented  in  references  12  to  15  and  19- 

For  some  cases  the  effects  of  aerodynamic  time  lag  are  important. 

There  are  two  different  sources  of  such  lag:  (1)  the  time  required  for 
an  aerodynamic  impulse  to  travel  from  one  component  of  the  airplane  to 
another  (for  example,  the  time  required  for  a  change  in  sldewash  at  the 
wing  to  reach  the  tail  -  a  phenomenon  commonly  referred  to  as  lag  of 
sidewash);  and  (2)  the  time  required  for  the  growth  and  decay  of  the 
aerodynamic  loads  on  the  airplane  components.  For  both  of  these  cases 
the  time-lag  effects  usually  become  increasingly  important  as  the  period 
of  the  lateral  oscillation  decreases.  The  effects  of  the  first  type  of 
time  lag  can  be  accounted  for  in  some  cases  by  modification  of  the  sta¬ 
bility  derivatives.  For  example,  the  effect  of  the  lag  of  sidewash  on 
the  derivative  Cn^  is  discussed  subsequently  under  the  section  on 

"Estimation  of  Lateral  Stability  Derivatives".  In  many  cases,  how¬ 
ever,  both  types  of  time  lag  will  require  special  stability  equations. 

No  general  treatment  of  these  cases  has  been  published  but  an  indi¬ 
cation  of  the  method  of  treatment  may  be  obtained  from  the  treatments 
of  autopilot  lag  in  references  13  nnd  15«  ^ 


CALCULATION  OF  PERIOD  AND  DAMPING 


As  pointed  out  in  references  1  and  2,  the  period  and  damping  of  the 
various  modes  of  the  lateral  motion  may  be  calculated  from  the  roots  of 
the  characteristic  equation 

AX^  +  BA,3  -h  CX^  +  DX,  +  E  =  0 


hy  the  equations 


P  = 


2jt 

I 


T 


and 


logp  2 


1/2 


0.693 
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where  R  represents  a  real  root  X  or  the  real  part  of  a  complex  root 
X  =  R  ±  li  and  I  represents  the  imaginary  part  of  a  complex  root. 
Negative  values  of  "^±^2  ^^P^^sent  the  time  required  to  double  amplitude 

for  unstable  modes  of  the  motion. 

The  values  of  the  coefficients  C,  D,  and  E  may  be 

obtained  by  the  method  given  in  steps  1,  2^  and  3  of  the  section  on 
"Calculation  of  Motions".  If  the  period  and  time  to  damp  are  to  be  cal¬ 
culated  for  a  number  of  related  cases^  however,  the  values  of  the  coef¬ 
ficients  A,  B,  C,  D,  and  E  may  be  more  conveniently  calculated  by 
a  tabular  procedure  such  as  that  shown  as  table  I  for  making  boundary 
calculations. 

Methods  of  determining  the  roots  of  the  biquadratic  characteristic 
equation  are  presented  in  appendix  C. 


CALCULATION  OF  MOTIONS 


Calculation  of  the  lateral  motions  of  an  airplane  involves  the 
integration  of  three  simultaneous  differential  equations  (see 
appendix  A)  to  obtain  a  general  solution  in  terms  of  the  mass  and 
aerodynamic  parameters  of  the  airplane.  The  general  equations,  once 
obtained,  can  then  be  used  to  obtain  numerically  the  motions  of  any 
airplane  in  terms  of  the  variation  with  time  of  the  angles  of  bank,  yaw, 

*  and  sideslip  or  some  function  of  these  angles  such  as  rolling  or  yawing 
velocity.  Various  methods,  such  as  those  given  in  references  20  to  22, 
are  of  course  available  for  integrating  the  differential  equations. 

w  Since  the  problems  met  in  airplane  dynamics  are  fairly  complex,  however, 

many  of  these  methods  are  not  suitable  because  of  the  difficulties  of 
computation  that  arise.  The  method  given  in  reference  4  (based  on  the 
Heaviside  operational  calculus)  is  satisfactory  for  calculating  the 
forced  motions  following  application  of  external  forces  or  moments  but, 
without  modification,  this  method  cannot  be  used  to  calculate  the  motions 
resulting  from  initial  displacements  in  bank,  yaw,  or  sideslip  or  from 
initial  values  of  rolling  6r  yawing  angular  velocity.  A  solution  based 
on  the  Laplace  transformation  is  more  satisfactory  than  that  based  on  the 
Heaviside  operational  calculus  because  it  permits  direct  calculation  of 
the  free  motions  following  any  initial  condition,  in  addition  to  calcu¬ 
lation  of  the  forced  motions  following  application  of  external  forces  and 
moments.  The  application  of  the  Laplace  transformation  to  the  calculation 
of  lateral  motions  is  outlined  in  appendix  B.  The  material  presented  in 
this  appendix  is  similar  to  the  work  presented  in  references  5  and  6 
except  that  the  mass  and  aerodynamic  stability  derivatives  have  been  com- 

*  bined  as  shown  in  appendix  A  to  reduce  the  number  of  arithmetical  and 
algebraic  processes  required  in  numerical  solutions. 
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The  process  of  calculating  the  motions  is  presented  as  a  series  of 
simple  though  lengthy  arithmetical  and  algebraic  steps  so  that  an  under¬ 
standing  of  the  calculus  involved  in  solving  the  differential  equations  ^ 

is  not  required.  The  method  as  shown  is  suitable  for  calculating  the 
motions  as  variations  of  0,  |3,  p,  and  r  with  time  for  the  case 

of  the  free  motions  following  initial  angular  displacements  (0o, 
and  Pq)  and  angular  velocities  (D0)q  and  (rii)o  case  of 

the  forced  motions  resulting  from  constant  impressed  forces  and  moments 
(Lc,  Nc,  and  Y,-).  These  are  the  cases  for  which  motions  are  usually 
calculated.  It  is  also  possible  to  calculate  the  motions  resulting  from 
impressed  forces  and  moments  which  are  arbitrary  functions  of  time  by 
the  methods  explained  in  references  6  and  7. 


Motions  Resulting  from  Initial  Angular  Displacements  and  Angular 

Velocities  and  from  Constant  Impressed  Forces  and  Moments 

The  six  steps  involved  in  obtaining  a  specific  solution  for  the 
lateral  motions  of  an  airplane  are: 

Step  1:  Determine  values  of  the  following  parameters: 

(a)  Mass  characteristics: 

^Zq,  ti,  and  p 

(b)  Geometric  characteristics: 

S  and  b 

(c)  Flight  conditions: 


V,  Cl,  and  7 

(d)  Aerodynamic  stability  derivatives: 


Cr,  ,  and 
■r 


The  methods  of  determining  the  values  of  the  aerodynamic  stability 
derivatives  are  given  in  subsequent  sections  of  this  paper. 


In  cases  where  impressed  forces  and  moments  are  used  as  disturbances, 
determine  the  values  of  the  factors 


that  are  appropriate  to  the  particular  problem. 
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y 


step  2:  From  the  known  factors,  evaluate  the  following  parameters 
which  are  the  stability  derivatives  in  the  form  in  which  they  are  used 
in  the  calculation  of  motions: 


K. 


I 


P 


I 


r 


2Kx2 


4k/ 


1 


K  - 

^2  IT  2 

m 

“  pSV 

C  T 

h 

■  2Kz2 

yp  = 

CZp 

■  i,K/ 

yp  = 

"  4K72  S 

Yr  = 

2  % 
4  % 


m 

pSb 


Also,  when  impressed  forces  and  moments  are  used,  evaluate 


= 


2K, 


2"Zc 


■n  —  P 

^  "  2k/  ^ 


yc  =  I  ^Y, 


p  2 

The  values  of  ^  and  determined  from  the  following 

expressions 


where 


2  2  2  2  2 
Kx  =  KXq  cos  t]  +  KZq  sin 

=  Kx^^cos^tj  +  Kx^^sin^T] 

KXZ  =  -  %o^)  ^ 


kZo 

b 


Step  3:  Solve  for  the  values  of  the 
following  coefficients  from  equations  (l) 


appropriate  ones  of  the 
to  (4): 


w 


llj- 
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In  all  cases  solve  for  the  values  of  C,  and  E 


where 


A  =  1  -  K1K2 

B  =  Pi  -  Ayp 

c  =  -Piyp  +  ^2  +  P5yp  +  -  ^6 


Cl  Cl  ^ 

D  =  P5  -^  +  P6  —  tan  7  +  Py 


Cl  Cl  ^ 

E  =  P3  -^  +  P4  tan  7 


1*1  =  -^p  "  “r  +  ^l“p 


Pg  =  ZpUj.  -  Z^Hp 


P3  =  Zpn^  -  Zj.np 


Pli  =  ^p^p  -  ^p°p 


P5  =  Kinp  -  Zp 


P6  =  KgZp  -  np 


P7  =  -P2yp  +  Psyp  +  J’i^yr  -  ^4 

The  quantities  Pi  to  F-j  are  factors  of  the  coefficients  B,  C,  D, 

and  E  which  are  comhinations  of  terms  that  occur  frequently  in  calcu¬ 
lations  of  motions  resulting  from  initial  angular  displacements  and 
velocities  and  which  are  consequently  grouped  together  for  convenience. 
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Calculate  the  values  of  Uq,  a^^,  .  .  .  a^  when  solving  for  the 
angle  of  bank  0  or  the  rolling  velocity  p: 

ao  =  ^^o^  1 


ai  = 

a^  =  -^qP^  +  (D0)o^-Ayp  +  KqI^  -  n^j  - 

(D>|f)o(Kinr  -  +  ^c  ’  ^c^l 


®-3 


8-4 


!^o(P6  ^  tan  7  +  Py)  -  ^  tan  7  -  +  (00)^  (  Pgy^. 

^6  -  K2Zryp  +  %yp)  +  (^)o(-Vr  "  V3)  " 

^c(“r  +  yp)  +  ^^(Kiyp  +  ^r)  -  ycP5 

(00^4  -  ^0^3  -  (Dlr)oP5)|i  tan  7  + 

^c(^P  -  +  "r^f)  /  ^c(^3yr  '  '  ^r^p)  '  ^c^S 


(2) 


a. 


=  ( 


tan  7 


y 


i6 
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Calculate  the  values  of  bQ,  b^,  •  •  ^5  when  solving  for 

angle  of  yaw  ''|f  or  the  yawing  velocity  r: 

bo  =  ^l^oA  ^ 

bi  =  >lroB  +  (Dt)oA 

”2  =  V  ■  '^o"’6  -  W'ol'Sj'p  -  “p)  ♦  *  ''l”p  -  ■ 

ZcK2  +  n^. 

b3  =  -(ZS^Pg  ^  +  Mro(p5  I"  +  P7)  -  P0P4  +  (D0)o(-P6yp  +  " 

Hpyp)  +  (°'^)o(Vp  ■  WP  ^  'p^p)^  'c(Vp  “p)  -  ^ 

'^c(^p  +  yp)  • 

t),  =  +  V3  -  (l^i^)oP6  +  +  ^c('ipyp  - 

0 

'^yp)  +  ^c  (^pyp  ■  ^p^i)  ■ 

=  (^c“p  ■  “c^p)t 
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Calculate  the  values  of  Cq, 
angle  of  sideslip  p: 


,  vhen  solving  for  the 


Cq  =  PqA 


Cl  =  0qA  +  tgA  tan  7  +  PqP^  +  (D0)oAy  -  (D>lf)QA(yj.  -  l)  +  y^A 


<=2  =  S^o^l  T  ^  Vl  +  Po^2  P  T  -  ^  ^ 


^1"P  -  Vr 


Vr  ■  ''p  +  (^2^r  -  •"  tan  7  +  Kinpy^  - 

Kl”p  -  Vr  ■"  *p  -  (’'l"r  -  *  'c  ("Vr  "■  '^2  *  + 

“c(yr  -  1  +  ycPl 

Cl  Cl  ,  /  Cl  Cl 

C3  =  00^2  “  +  ^0^2  —  tan  7  +  (D0)o( -KgZp  tan  7  +  np  ^  tan  7 


K2^r  ”  “  y  +  (l>^)o{Kinp  ~  tan  7  -  Zp  ^  tan  7 


^1%"  2  ^  2/  ^  ”  ^^p  ^2  ~  ^2  2  ^  ^ 

^  Cl  Cl  \ 

^ctVr  ^  ^ryp  -  %  —  +  —  tan  ^  +  y^Pg 


tan  rj 


%  =  ^  tan  7  -  rv  (Z^  -  Z^ 


Step  4:  Solve  for  the  roots  \2>  ^3^  ^4 

biquadratic  equation 

AX^  +  Bx3  +  CX^  +  dx  +  e  =  0 
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where  the  values  of  the  coefficients  A,  ^  etc.  were  given  by 

the  solution  of  equations  (l).  Methods  of  determining  the  roots  of  the 
biquadratic  equation  are  given  in  appendix  C. 

Step  5:  Use  the  coefficients  obtained  from  equations  (l)  to  (4) 
and  the  roots  of  equation  (5)  to  solve  for  the  following  coefficients: 

Calculate  the  values  of  the  factors  A]_,  A2^  .  .  .  ^  Ag  when 

solving  for  the  angle  of  bank  0  or  the  rolling  velocity  p: 


Ai  = 


^0^1^  ^1^1^  ^*2^1^  ^*3^1^  ®*4^1 

SAX.^  +  5BX.i^  +  +  2EXi 


Ao  = 


6aX2^  +  5BX2^  +  *tCX2^  +  SDXg^  +  2EX2 


A3  = 


6AX^  +  56X3^  +  40X3^  +  30X3^  +  2EX3 


^■0^4^  +  a2^X]|_^  +  a2X||_^  +  a3X)|_^  +  a]|_X]|_  +  a^ 
^4  ”  R  li  P 

6aX1(.^  +  5BXi^  +  4CX4^  +  3DX4  +  2EX4 


A 

^5  E 


^6  =  i-K  - 


D 

^5  E 
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Calculate  the  values  of  the  factors  B2,  .  .  Bg  when 

solving  for  the  angle  of  yaw  or  the  yawing  velocity  r: 


Bi  = 


_  Bpli^  +  ^1^1^  +  hgX]^^  + 


6aX.i5  +  5BX,2^^  +  itCXj^S  +  3dx,^2  + 


B  =  ^0^2^  +  ^1^2^  •*•  ^2^2^  +  ^3^2^  ^4^2  + 

6AX,2^  +  5BX2^  +  ^^2^  +  30X2^  +  2EX2 


B3  = 


^0^3^  ■*■  ^1^3^  ■*■  1*2X2^  +  B^X^^  +  bj^X^  +  b^ 


6aX3^  +  5BX3^  +  40X3^  +  3BX3^  +  2EX3 


Bi,  = 


^0^4^  +  bj^Xij^^  +  b2X]j^^  +  b3X|j^^  +  b]j^X)|  +  b^ 

SAXlp  +  5BX4^  +  4cX4^  +  3DX4^  +  2EX4 


i 


« 


(7) 
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Calculate  the  values  of  the  factors  Cg, 

solving  for  the  angle  of  sideslip  P: 


.when 


CqX.!^  +  +  <=1<-^1 

6aXi^  +  5BX3_^  +  4CXi^  +  3DXj^^  +  2EXi 

CpXg^  +  c^x/  +  CgXg3 

6AX2^  +  5BX2^  +  ^^-2^  +  BDXg^  +  2EX2 


"3 


C0X35  +  C^Xg^  +  CgXc^S  +  C^Xg^  +  C4X3 

6AX^^  +  55X3^  +  ^X33  +  3DX3^  +  2EX3 


(8) 


CqX^^  +  <^2^X4^  +  02X^3  +  Ci^X^^  +  cj^Xij. 

6ax4^  +  5BX4^  +  40X4^  +  30X4^  +  2EX4 


"5  J 

If  equation  (5)  has  conjugate  complex  roots,  the  values  of  the 
coefficients  (equations  (6)  to  (8))  corresponding  to  these  roots  'wiH 
be  conjugate  complex.  In  order  to  facilitate  treatment  of  this  case  it 
is  convenient  to  establish  some  special  notation.  This  special  notation 
is  explained  in  appendix  D. 


I 


if 
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Step  6:  The  equations  of  motion  are  written  in  different  form 
depending  upon  the  roots  of  equation  (5).  If  the  characteristic  equa¬ 
tion  has  four  real  roots  “the  general  form  of 

the  equations  of  motion  is  used,  as  follows: 

I  crX"]  cXn  ctXo 

0  =  Aie  ^  ■+  Age  ^  +  A3e  ^  +  Aj^e  ^  +  A^a  +  Ag 

crXi  oXn  oXo  crXji 

^  =  B^e  +  B2e  ^  -f  B^e  ^  +  B]^e  ^  +  B^a  +  Eg 

oX-\  0X0  0X0  0X)i  , 

3  =  C^e  -L  +  Cge  2  +  c^e  3  +  c^e  +  C5  >  (9) 

If  *  cX-i  0X0  0X0  crX),  \ 

p  =  ^Aj^X^^e  ^  ^ 

If  0X1  0X0  0X0  ^X|i  \ 

r  =  ^^1^1^  B^X^s  +  ^  ■*"  ^4^4^  ^ 

If,  as  is  generally  the  case,  equation  (5)  has  two  complex  roots  and  two 
real  roots  (R  +  li,  R  -  li,  X^,  and  X^^),  the  equations  of  motion  may  be 
expressed  as 

/  0R  ,  .  0X0  0X), 

p  =  K^e  cos(0l  +  cjD^)  +  A^e  -j  +  A]^e  ^  +  A^0  +  Ag 

cos  (0I  +  a>g)  +  B^e^^S  +  B^^e^^^  +  B^0  +  Bg 
p  =  cos  (0I  +  cOq)  + 

P  =  ^|ka|/r^  +  I^  cos/01  +  coa  +  tan^^ 


“  +  I^  e^^  cos^I  +  co^  +  tan^^ 

+  Ai^X^e^^^  +  A^ 

=  +  I^  e^^  cos^I  +  a>g  +  tan"^  ^ 

0X0  o^Xji  ”1 

B^X^G  ^  +  Bj^Xj^e  +  B^j 


4 
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where 


“>a' 

=  tan**^ 

Ia' 

1 

%' 

=  2I 

a>B' 

=  tan”^ 

Ib’ 

%’ 

(lla) 

Kc’ 

=  2\ 

4c +  Ic’2 

=  tan"^ 

Ic’ 

Rp' 

The  coefficients 

'  <^A> 

and 

cjoq  are  defined 

in  equa- 

tions  (lOa)  and 

Ia,  Ra'^  ^a' 

are 

defined 

in  appendix  D. 

Solve  the  appropriate  ones  of  these  equations  of  motion  (equa¬ 
tions  (9)>  (10)^  or  (11))  by  substituting  values  of  the  nondimen sional 
time  factor  a  in  the  equations  and  solving  for  ^  ^  or  r. 


Motions  Resulting  from  Arbitrary  Disturbances 

The  motions  resulting  from  arbitrary  forcing  functions  can  be 
obtained  from  the  motions  resulting  from  constant  impressed  forces  and 
moments  by  the  methods  explained  in  references  6  and  7. 

A  very  useful  method  of  obtaining  the  motion  resulting  from  various 
abrupt  gust  and  control  disturbances  is  given  by  Jones  in  reference  7* 

In  this  paper  it  is  pointed  out  that,  although  the  component  motions  of 
an  airplane  must  be  calculated  simultaneously  (that  is,  by  simultaneous 
differential  equations),  the  effects  of  component  disturbances  may  by 
the  principle  of  superposition  be  calculated  separately  and  later  added 
in  any  desired  proportion.  Thus,  if  a  given  rolling  moment  causes  a 
20°  bank  in  1  second  and  if  a  given  yawing  moment  causes  a  5^  bank  in 
1  second,  the  combined  effect  of  both  acting  simultaneously  will  be  a 
25°  bank  in  1  second,  Jones  also  points  out  a  somewhat  similar  fact 
with  regard  to  the  effects  of  disturbances  that  are  not  applied  simul¬ 
taneously.  This  fact  is  that,  if  a  given  disturbance  which  arises  at 
the  time  t  =  0  is  later  augmented,  the  effect  of  the  increment  of 
disturbance  will  run  its  course  independently  of  the  effect  of  the 
original  disturbaiace ,  For  example,  in  a  problem  involving  the  correc¬ 
tion  for  a  gust  disturbance  by  a  manipulation  of  the  control,  the  motion 
produced  by  the  gust  disturbance  can  be  calculated  independently  and  the 
motion  caused  by  the  assumed  corrective  control  manipulation  can  be 
added  to  it  at  any  desired  point.  This  example  is  illustrated  graphi¬ 
cally  in  figure  3. 
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The  principle  of  superposition  may  he  applied  analytically  as  well 
as  graphically.  The  analytical  application  which  makes  use  of  Carson's 
integral  or  Duhamel's  integral  is  described  in  references  7  23. 

This  method  is  useful  for  calculating  the  motions  resulting  from 
impressed  forces  and  moments  which  are  arbitrary  functions  of  time.  By 
application  of  these  methods^  the  solutions  for  constant  impressed 
forces  and  moments  can  be  used  to  obtain  new  solutions  for  any  arbitrary 
variation  of  impressed  forces  and  moments  with  time  which  can  be 
expressed  by  a  mathematical  formula.  Some  simple  variations  of 
impressed  forces  and  moments  with  time  and  their  Laplace  transforms  are 
given  in  reference  6.  The  transforms  for  any  other  function  for  which 
transforms  have  been  worked  out  may  be  found  in  tables  of  Laplace 
transforms . 


CALCULATION  OF  STABILITY  BOUNBARIHIS 
Oscillatory  Stability  Boundaries 


As  pointed  out  in  the  preceding  section  of  this  report^  the  degree 
of  stability  of  the  uncontrolled  motions  of  an  airplane  is  indicated  by 
roots  of  the  characteristic  equation 

AX^  +  BX^  +  CX^  +  DX  +  E  =  0 

For  stability  the  real  roots  or  the  real  part  of  the  complex  roots  of 
the  characteristic  equation  must  be  negative.  A  useful  discriminant 
for  determining  some  of  the  characteristics  of  the  roots  in  stability 
work  is  Routh's  discriminant  R  (R  =  BCD  -  AD^  -  B^E),  The  use  of  this 
discriminant  in  dynamic  stability  analyses  has  been  pointed  out  in  many 
reports^  for  example^  references  1,  2^  3^  5^  21^  and  24.  Routh  has 
shown  (reference  20)  that,  if  R  and  the  coefficient  E  are  finite^ 
the  necessary  and  sufficient  conditions  that  the  real  roots  and  the 
real  parts  of  the  complex  roots  should  be  negative  are  that  every  coef¬ 
ficient  of  the  biquadratic  and  also  R  should  have  the  same  sign. 

Routh  also  showed  that  when  R  =  0  and  B  and  D  have  the  same  sign 
there  are  a  pair  of  complex  roots  with  the  real  parts  zero.  Since  the 
value  of  the  real  part  of  a  complex  root  indicates  the  stability  of  an 
oscillatory  mode  of  the  motion  of  an  airplane,  the  lateral  oscillation 
is  neutrally  stable  when  R  =  0  and  the  coefficients  B  and  D  have 
the  same  sign.  Oscillatory  stability  boundaries  can  be  determined, 
therefore,  by  solving  the  equation  R  =  0  and  checking  to  determine 
whether  the  signs  of  B  and  D  are  the  same. 

Since  two  of  the  most  important  stability  derivatives  affecting 
lateral  stability  are  the  directional  stability  derivative  Cn^  and 
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the  effective  dihedral  derivative  boundaries  for  neutral  oscil- 

latory  stability  are  usually  calculated  as  a  function  of  these  two 
derivatives  as  illustrated  in  figure  4.  These  calculations  are  gener¬ 
ally  carried  out  by  the  method  shown  in  table  I.  This  table  contains  a 
numerical  example  and  step-by-step  instructions  for  using  the  table. 

The  results  of  this  numerical  example  are  plotted  in  figure  k.  The 
procedure  illustrated  in  table  I  is  first  to  assume  values  of  the  inde¬ 
pendent  variable  Cnp  to  cover  the  range  for  which  the  boundary  is 
required.  The  values  of  all  the  other  mass  and  aerodynamic  stability 
derivatives  except  Czp  are  then  estimated.  The  value  of  Cnp  is 
generally  assumed  to  have  been  varied  by  varying  the  size  of  the  verti¬ 
cal  tail  and  consequently  the  tail  contribution  to  each  of  the  other 
stability  derivatives  varies  as  Cnp  is  varied.  The  values  of  the 
coefficients  and  E  and  then  R  are  calculated  as 

functions  of  Zp! 


The  values  of  Zp  corresponding  to  the  assumed  values  of  Cnp  for  the 

condition  of  neutral  oscillatory  stability  are  next  obtained  by  solving 
the  expression  R  =  0  which  is  a  quadratic  in  Zp  that  is  of  the  form 

p 

^1^3  '^l^P  +  w^  =  0 

Finally^  the  values  of  C^^  corresponding  to  the  assumed  values  of  Cnp 
are  obtained  from  the  values  of  Zp. 

The  values  of  Zp  which  satisfy  the  expression  R  =  0  must  be 
checked  to  determine  whether  they  satisfy  the  other  condition  for 
neutral  oscillatory  stability  -  that  the  sign  of  the  coefficients  B 
and  D  must  be  the  same.  This  check  can  be  performed  readily  by  sub¬ 
stituting  the  values  of  Zp  which  satisfy  R  =  0  into  the  expression 
for  D  which  is  a  linear  equation  of  the  form 


D  =  u2Zp  +  Vp 


Thus,  the  sign  of  D  is  determined.  The  sign  of  B  is  a  constant  for 
any  given  value  of  Cnp  and  is  almost  invariably  positive  since  the 


three  predominant  terms  of  B  contain  the  derivatives  C?  , 


and  Cyp  which  in  all  practical  cases  contribute  a  positive  increment 
to  the  value  of  B. 
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Since  t-vo  values  of  Czp  satisfy  the  condition  R  =  0  for  each 
value  of  Cup,  the  R  =  0  curve  has  two  branches.  As  pointed  out  in 

reference  24;,  one  of  the  branches  of  the  R  =  0  curve  generally  repre¬ 
sents  an  oscillatory  stability  boundary  and  the  other  branch  represents 
a  line  of  numerically  equal  real  roots  with  opposite  signs.  (See 
fig.  4.)  If  neither  of  the  values  of  which  satisfy  the  expres¬ 

sion  R  =  0  for  a  particular  value  of  Cnp  is  found  to  represent  a 

point  of  neutral  oscillatory  stability,  the  lateral  motion  has  no  oscil¬ 
latory  mode  for  that  value  of  0^^.  If  l^oth  of  the  values  of 

which  satisfy  the  expression  R  =  0  are  found  to  represent  points  of 
neutral  oscillatory  stability,  the  lateral  motion  has  two  oscillatory 
modes.  In  this  case,  since  the  boundary  D  =  0  represents  the  line  of 
infinite  period,  the  branch  of  the  R  =  0  boundary  which  lies  close  to 
the  D  =  0  boundary  is  usually  the  boundary  for  neutral  stability  of 
the  longer  period  of  the  two  oscillatory  modes.  A  detailed  discussion 
of  the  significance  of  the  stability  boundaries  and  the  regions  formed 
by  these  boundaries  is  given  in  reference  24. 

In  calculating  stability  boundaries  for  a  specific  airplane  a  com¬ 
plete  solution  such  as  that  explained  in  the  preceding  paragraphs  should 
be  made.  For  general  studies  of  stability,  however,  approximate  oscil¬ 
latory  stability  boundaries  may  be  calculated  much  more  simply  by  the 
methods  shown  in  reference  24. 

As  pointed  out  previously,  methods  of  calculating  lines  of  constant 
period  and  damping  of  the  lateral  oscillation  are  presented  in  refer¬ 
ences  8  and  9" 


Spiral  Stability  Boundaries 

Spiral  stability  boundaries,  like  oscillatory  stability  boundaries, 
are  usually  determined  as  a  function  of  the  directional  stability  deriva¬ 
tive  Cnp  and  the  effective  dihedral  derivative  as  illustrated  in 

figure  4.  As  pointed  out  in  reference  1,  neutral  spiral  stability 
occurs  when  the  E  coefficient  of  the  characteristic  equation  is  zero 
(E  =  0).  A  spiral  stability  boundary  can  be  easily  obtained  from  this 
relation.  If  expressions  for  E  (in  terms  of  Zp)  corresponding  to 
several  values  of  Cnp  have  already  been  obtained  in  the  process  of 

calculating  an  oscillatory  stability  boundary,  the  equations  formed  by 
setting  these  expressions  for  E  equal  to  zero  can  be  solved  for  'the 
values  of  Zp  ^and  hence  C^p)  corresponding  to  the  assumed  values 

of  C„  .  If  the  values  of  E  have  not  already  been  obtained  in  the 
process  of  calculating  an  oscillatory  stability  boundary,  a  spiral 
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stability  boundary  for  the  level-flight  condition  (7  =  0)  can  be  cal¬ 
culated  simply  from  the  equation 


Clr 

-  Cn^  Cnp 


(12) 


Values  of  Cnp  are  assumed  within  the  range  for  which  the  boundary  is 
required.  The  values  of  and  corresponding  to  each  value  of 

Cnj3  then  determined.  The  tall  contributions  to  these  derivatives 

generally  vary  with  Cnp  since  Cn^  is  usually  assigned  to  be  varied 
by  changing  the  size  of  the  vertical  tall. 


ESTIMATION  OF  LATERAL  STABILITY 

DERIVATIVES 


GENERAL  REMARKS 


Methods  of  estimating  the  lateral  stability  derivatives  have  been 
presented  in  numerous  publications  but  no  single  report  has  contained 
information  for  estimating  the  contribution  of  all  principal  airplane 
components  to  all  the  derivatives  for  airplanes  having  ajiy  sweep  angle 
or  aspect  ratio.  In  the  present  paper,  an  approach  to  such  a  presenta¬ 
tion  is  made  by  the  coordination  of  and  reference  to  existing  estima¬ 
tion  methods,  by  reference  to  publications  containing  data  which  should 
be  useful  in  making  estimates,  and  by  the  suggestion  in  some  cases  of 
simple  new  empirical  formulas.  Detailed  estimation  methods  are  pre¬ 
sented  for  low-subsonic -speed  conditions  but  only  a  brief  discussion  and 
a  list  of  references  are  given  for  transonic-  and  supersonic -speed  con¬ 
ditions.  In  general,  the  estimation  methods  presented  should  be  expected 
to  yield  only  fairly  accurate  values  suitable  for  making  first  approxi¬ 
mations  of  dynamic  stability.  This  limitation  applies  especially  to 
the  cases  in  which  the  derivatives  are  based  completely  on  theory. 

For  convenience,  the  references  that  should  be  useful  in  estimating 
the  stability  derivatives  are  presented  in  table  II.  The  references 
are  grouped  according  to  the  speed  range  covered  ( subsonic  or  super¬ 
sonic)  and  according  to  the  derivatives  presented  in  each  report.  The 
references  for  the  subsonic  case  (references  1  and  25  to  94)  are  further 
divided  into  two  groups  -  one  including  reports  which  contain  estimation 
methods  and  the  other  including  reports  which  contain  experimental  data 
that  should  be  useful  in  making  estimates  of  derivatives.  The 
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references  for  the  supersonic  case  (references  95  "to  115)  ar®  sub¬ 
divided  according  to  wing  plan  form. 

The  following  sections  covering  the  estimation  of  the  nine  sta¬ 
bility  derivatives  are  divided  into  three  groups  according  to  the  tpe 
of  derivative  -  sideslip  derivatives  ^np^  deriva¬ 

tives  (Cnp.  C^p,  Cypj,  and  yawing  derivatives  (Cnr,  Clr>  ^Yr)  • 
derivatives  Gyp  and  Cy^.  have  usually  been  neglected  in  making 

dynamic  lateral  stability  calculations  because  theory  indicated  that  for 
unswept  wings  Gyp  and  Gy^.  were  zero.  Recent  experimental  data, 

however,  have  indicated  that  both  swept  and  unswept  wings  produce  meas¬ 
urable  values  of  these  derivatives  (references  25,  59,  and  86).^  Since 
the  vertical  tail  contributes  to  Gyp  and  Gy^.,  it  appears  desirable 

to  estimate  these  derivatives  and  to  use  them  in  the  calculations  of 
stability  unless  it  is  established  that  for  the  case  in  question  the 
effects  of  Gyp  and  Gy^.  on  stability  are  negligible.  For  these  two 

derivatives,  only  the  effect  of  the  wing  and  vertical  tail  need  to  be 
considered. 

The  methods  of  estimating  the  rolling  and  yawing  derivatives  pre¬ 
sented  herein  were  obtained  from  theoretical  treatments  based  on  the 
assumption  of  steady  rolling  and  yawing  and  from  experimental  data  _  . 

obtained  principally  from  tests  made  under  conditions  of  steady  rolling 
and  yawing.  The  only  information  that  applies  directly  to  the  oscil-^ 
latory  case  is  a  limited  amount  of  data  on  Gn^  obtained  by  oscillation 

techniques.  When  calculations  are  made  in  which  the  oscillatory  mode 
is  the  subject  of  interest,  some  consideration  should  be  given  to  cor¬ 
recting  the  derivatives  based  on  steady  rolling  or  yawing  to  account 
for  differences  in  the  derivatives  that  are  likely  to  exist  as  a  result 
of  differences  between  the  oscillatory  motion  and  the  steady  rolling 
and  yawing  motion.  For  example,  the  data  of  reference  82  have  indi¬ 
cated  that,  for  flap-extended  or  power-on  conditions,  fairly  large  dif¬ 
ferences  mi^t  exist  between  the  values  of  the  tall  contribution  to  Gnj- 
for  the  steady  yawing  and  yawing  oscillation  cases.  At  present  little 
information  is  available  for  correcting  the  values  of  Gn^.  for  fhe 

steady  yawing  case  to  apply  to  the  oscillatory  case  and,  unfortunately, 
little  or  no  information  is  available  for  correcting  the  other  stability 
derivatives. 

Since  most  wind-tunnel  force-test  data  that  are  likely  to  be  used 
in  making  estimates  of  the  stability  derivatives  are  probably  for  much 
lower  Reynolds  numbers  than  those  for  the  full-scale  airplane,  some 
adjustments  to  the  data  are  usually  required  to  account  for  the  dif¬ 
ferences  in  Reynolds  number.  The  effects  of  Reynolds  number  should  be 
considered  in  the  cases  of  all  the  derivatives,  especially  those  which 
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are  estimated  by  methods  that  involve  the  use  of  force-test  data. 

Methods  of  correcting  for  Reynolds  number  effects  for  some  of  the 
derivatives  are  discussed  in  the  following  sections  which  cover  the 
estimation  procedures.  In  the  cases  where  the  Reynolds  number  effects 
are  not  discussed^  it  can  be  assumed  that  any  abrupt  variation  in  the 
derivatives  near  the  stall  for  low-scale  data  will  also  be  present  for 
the  full-scale  airplane  but  will  probably  occur  at  a  higher  lift  coef¬ 
ficient  because  of  the  higher  maximum  lift  coefficient  of  the  airplane. 
An  indication  of  the  lift-coefficient  range  over  which  the  theory  may 
not  be  expected  to  give  reliable  values  of  stability  derivatives  for  the 
full-scale  airplane  can  be  obtained  from  large-scale  drag  data.  The 
analysis  of  reference  86  indicates  that  the  variation  of  the  derivatives 
with  lift  coefficient  is  different  from  the  theoretical  variation  at 
lift  coefficients  above  that  at  which  the  drag  due  to  lift  increases 
abruptly  from  the  ideal  value  Cl^/ttA. 

The  effects  of  Mach  number  and  power  are  not  treated  in  the  sections 
on  the  individual  derivatives  but  are  discussed  briefly  in  separate 
sections.  A  detailed  treatment  of  these  effects,  including  design 
formulas  and  charts,  was  considered  beyond  the  scope  of  this  paper. 


THE  SIDESLIP  DERIVATIVES  CYo^  Cji  ,  Ci 

H  M  H 


No  satisfactory  purely  theoretical  methods  have  yet  been  developed 
for  obtaining  accurate  estimates  of  the  sideslip  derivatives  Cy^,  Cup, 

and  C^p  for  a  complete  airplane,  primarily  because  of  large  inter¬ 


ference  effects  between  the  various  airplane  components  and  because  of 
large,  and  often  unpredictable,  variations  of  the  derivatives  with  angle 
of  attack.  Fortunately,  these  derivatives  can  be  obtained  from  conven¬ 
tional  wind-tunnel  force -test  data.  Such  experimental  data  are  essential 
to  the  accurate  determination  of  sideslip  derivatives.  It  is,  of  course, 
highly  desirable  to  have  force-test  data  for  the  exact  airplane  design 
under  consideration,  but  reasonably  accurate  estimates  can  usually  be 
made  by  correcting  the  force-test  data  for  a  generally  similar  design. 

The  methods  of  correcting  the  force-test  data  on  a  similar  design  for 
use  in  the  case  under  consideration  are  covered  in  the  following  sec¬ 
tions.  In  the  formulas  presented,  the  subscript  word  "design”  is  used 
to  designate  the  design  under  consideration  and  the  subscript  word  "data" 
is  used  to  designate  the  similar  design  for  which  force-test  data  are 
available . 


Force-test  data  should  be  used  to  determine  the  effect  on  the  side¬ 
slip  derivatives  of  such  airplane  components  as  leading-edge  high-lift 
devices,  stall-control  devices,  trailing-edge  flaps,  nacelles,  external 
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stores,  canopies,  and  dorsal  and  ventral  fins.  The  effect  of  leading- 
edge  high-lift  devices  is  usually  merely  to  extend  to  a  higher  lift 
coefficient  the  same  variation  of  the  derivative  with  lift  coefficient 
as  for  the  plain  wing.  Trailing-edge  flaps  often  have  large  effects  on 
the  contributions  of  both  the  wing  and  the  vertical  tail  to  the  sideslip 
derivatives  (references  39  snd  69))  and  since  these  effects  are  not 
easily  estimated,  it  appears  that  in  these  cases  use  of  force-test  data 
is  essential.  The  addition  of  nacelles  and  external  stores  generally 
has  been  found  to  decrease  the  directional  stability  factor  Cnp 

slightly.  The  results  of  a  limited  amount  of  research  to  determine  the 
effect  on  the  sideslip  derivatives  of  the  size  and  shape  of  canopies 
has  been  reported  in  references  48  and  Y3  but  these  results  axe  inade¬ 
quate  for  making  accurate  predictions  of  the  effects  of  canopies.  The 
effects  on  the  sideslip  derivatives  of  dorsal  and  ventral  fins  are 
usually  small  at  the  small  and  moderate  angles  of  yaw  that  are  generally 
considered  in  stability  calculations.  (See  references  4-7  and  71- ) 


In  estimates  of  the  lateral  force  due  to  sideslip  derivative 

force -test  data  for  the  design  under  consideration  should  he  used  when¬ 
ever  possible.  If  such  data  are  not  available^  data  for  a  similar 
design  can  be  used  and  corrected  as  follows: 

Wing-fuselage . -  Since  the  wing-fuselage  contribution  to  Gyp  is 

usually  relatively  small  compared  with  that  of  the  vertical  tail^  great 
accuracy  is  not  rec^uired  in  estimating  this  factor.  This  contribution 
may  be  estimated  as  follows: 

(1)  Wing:  If  the  wings  of  the  two  designs  are  generally  similar 

the  difference  in  Cyn  can  be  considered  negligible  and  no  correc- 

Pwing 

tion  is  necessary.  The  theory  of  reference  25  does  not  appear  to  be 
suitable  for  use  in  estimating 

(2)  Fuselage:  If  the  two  fuselages  are  similar  in  shape^  the 

difference  in  probably  be  estimated  satisfactorily  by  cor¬ 

recting  for  the  difference  in  the  relative  size  of  the  fuselage  and 
wing  for  the  two  airplanes.  It  appears,  however,  from  table  X  of  refer¬ 
ence  69  unlikely  that  a  reliable  prediction  of  Cyp^^^  can  be  made 

directly  from  the  geometry  of  the  fuselage.  Some  additional  data  on 
CYq  nre  presented  in  reference  77*  Experimental  data  from  other 

Pfus 

investigations  have  shown  that  differences  in  fuselage  cross-section 
can  cause  very  large  differences  in  the  variation  of  Cyp^^^  with 

angle  of  attack.  For  example,  in  the  case  of  a  flat  fuselage  with  the 

major  cross-sectional  axis  horizontal,  the  sign  of  Cyn  has  been 

^fus 
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found  to  reverse  at  moderate  and  high  angles  of  attack.  Force -test 
data  are  essential  for  making  estimates  in  such  cases. 

(3)  Wing-fuselage  interference:  For  low-wing  or  high-wing  con¬ 
figurations^  wing-fuselage  interference  causes  the  value  of  Cyp  to  be 

greater  than  that  obtained  by  adding  the  contributions  of  the  wing  and 
fuselage.  (See  reference  39.)  If  the  vertical  location  of  the  wing  on 
the  fuselage  is  generally  similar  for  the  two  designs^  however^  any 
correction  for  a  difference  in  this  interference  factor  can  be  neglected. 

Vertical  tail.-  Accurate  estimates  of  necessary 

because  this  factor  is  used  to  estimate  the  tail  contribution  to  several 
other  derivatives.  This  factor  is  especially  important  at  low  angles 
of  attack  because  in  this  case  the  tail  contribution  is  often  much 
greater  than  the  wing-fuselage  contribution  to  all  derivatives  except 
C^p.  For  this  reason  it  is  highly  desirable  to  have  tail-off  and  tail- 
on  force-test  data  for  the  design  under  consideration  or  for  a  very 
similar  design.  Corrections  to  the  data  for  a  similar  design  can  be 
made  as  follows: 

(1)  Correction  for  differences  in  wing  area^  tail  area^  and  tail 
lift -curve  slope  can  be  made  by  the  following  formula: 


(^Yptail) 


design 


data 


Bdata 

^design 


(13) 


The  value  of  Cl^  can  be  obtained  from  figures  5  and  6  which  are 

ntail 

based  on  the  theory  of  reference  34  and  on  the  theory  and  data  of  refer¬ 
ences  28  and  35*  The  chart  of  figure  6  can  be  used  to  estimate  the 
change  in  the  effective  aspect  ratio  of  the  vertical  tail  caused  by  the 
end-plate  effect  of  the  horizontal  tail.  It  should  be  emphasized  that 
for  the  best  accuracy  the  charts  in  figures  5  and  6  should  be  used  in 
conjunction  with  formula  (13)  Tor  correcting  existing  force-test  data 
and  not  for  making  a  direct  estimate  of 


(2)  In  the  case  of  V-tails.  the  correction  for  Cvr  can  be 

Ptail 


made  as  follows: 


Pv-tail ' (iesi 


J^KC]^^Sv-tail  sin^rj 
Pv-tail)data  ^KCL^^Sy-tail  sin2rj 


design  data 

\  ^design 

/data  ^ 
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where  the  terms  CLa^ .  T,  and  K  are  the  same  as  given  In  reference  30 
and  are  defined  as  follows: 

Ct  slope  of  the  tail  lift  curve  in  pitch  measured  in  the  plane 

normal  to  the  chord  plane  of  each  tall  panel 

r  dihedral  angle  of  tail  surface  measured  from  XY -plane  of  the 

tail  to  each  tail  panel,  degrees 

K  ratio  of  sum  of  lifts  obtained  by  equal  and  opposite  changes 

in  angle  of  attack  of  two  semispans  of  tail  to  lifts  obtained 
by  an  equal  change  in  angle  of  attack  for  the  complete  tail 

Values  of  the  term  K,  which  are  usually  about  0.7^  can  be  obtained  from 
reference  30. 

(3)  Since  large  differences  in  sidewash  and  dynamic  pressure  at 
the  tall  can  be  caused  by  differences  in  wing  plan  form  and  wing  loca¬ 
tion,  use  of  experimental  data  for  the  specific  design  or  at  least  for 
a  design  which  has  a  closely  similar  wing-fuselage  combination  and 
vertical  tail  location  is  extremely  desirable.  No  methods  are  available 
which  permit  accurate  predictions  of  sidewash  at  the  tall,  but  the^ 
experimental  data  of  references  39)  ^9)  ^9  can  be  used  to  obtain 

some  indication  of  the  variation  in  sidewash  with  vertical  location  of 
an  unswept  wing  on  a  fuselage  and  the  experimental  data  of  references  3° 
and  77  provide  additional  Information  on  sidewash  at  the  tail.  Other 
experimental  data  indicate  that  the  sidewash  fields  produced  by  hi^ly- 
swept,  low-aspect -ratio  wings  or  by  fuselages  of  flat  cross  section  can 
sometimes  be  strong  enou^  at  high  angles  of  attack  to  reverse  the 
effectiveness  of  a  conventionally-located  vertical  tail  surface.  Until 
a  reliable  method  is  developed  for  predicting  these  large  sidewash 
effects,  force -test  data  appear  to  be  the  only  means  by  which  satisfac¬ 
tory  estimates  of  can  be  obtained. 


Cnp 

Although  attempts  have  been  made  to  develop  methods  for  estimattog 
the  yawing  moment  due  to  sideslip  (static  directional  stability)  deriva¬ 
tive  Cnp  (for  exanqile,  references  68  and  69)  no  reliable  method  has 

yet  been  obtained.  The  use  of  force-test  data  therefore  seems  imperative. 

Force-test  data  for  the  design  under  consideration  should  be  used 
if  available.  If  such  data  are  not  available,  use  data  for  a  similar 
design  and  correct  as  explained  in  the  sections  to  follow. 
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Wing-fuselage . ~  The  corrections  for  the  wing-fuselage  contributions 


(1)  Correction  f^  wing  -  From  figure  7  (taken  from  reference  25) 
the  values  of  ^Ciij3/CL  design  under  consideration  and  for 

the  design  for  which  test  data  are  available  can  be  determined.  The 
effect  of  differences  in  taper  ratio  can  be  neglected.  (See  refer¬ 
ences  6o  and  66.)  The  difference  between  these  values  of 

should  then  be  added  (with  proper  regard  for  sign)  to  the  experimental 
data  for  the  complete  model. 


(2)  Correction  for  fuselage  -  The  formula 


can  be  used  to  calculate  the  Cnp  of  the  fuselage  (per  radian)  for  the 

design  under  consideration  and  for  the  similar  design  for  which  force- 
test  data  are  available.  The  differences  between  these  two  values  can 
then  be  added  (with  proper  regard  for  sign)  to  the  force -test  data  for 
the  complete  model.  Formula  (15)  does  not  include  the  effect  of  fine¬ 
ness  ratio  and  should  not  be  used  for  fineness  ratios  less  than  4.  This 
formula  is  an  approximate  empirical  expression  which  should  not  be  used 

to  estimate  the  value  of  Cun  directly  but  should  only  be  used  as 

Pfus 

indicated  to  determine  a  correction  for  force-test  data.  This  correc¬ 
tion  method  should  not  be  used  in  the  cases  of  high  angles  of  attack 
when  there  are  large  differences  in  fuselage  configuration.  Force-test 
data  are  essential  in  such  cases. 


(3)  Correction  for  vertical  location  of  the  wing  -  If  the  designs 
are  generally  similar^  the  correction  for  the  vertical  location  of  the 
wing  on  the  fuselage  can  be  neglected.  (See  reference  39.) 

(k)  Correction  for  center -of -gravity  position  -  If  the  center-of- 
gravity  position  for  the  design  under  consideration  is  appreciably  dif¬ 
ferent  from  that  for  the  design  for  which  force-test  data  are  available, 
the  value  of  Cnp  for  the  wing-fuselage  combination  can  be  corrected 

by  multiplying  the  value  of  Cyp  for  the  wing-fuselage  combination  by 

the  distance  between  center -of -gravity  positions  (expressed  in  wing 
spans ) . 
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Vertical  tall.-  Corrections  to 


Cyptaii  length  z/h 


can  be  made  by 


for  differences  in 
the  following  formula: 


The  contribution  of  wing-tip  fins  to  Cnp  is  treated  in  refer¬ 
ences  70  84, 


In  estimates  of  the  rolling  moment  due  to  sideslip  (effective 
dihedral)  force -test  data  for  the  design  under  consideration  should 

be  used.  If  such  data  are  not  available,  data  for  a  similar  design  can 
be  used  and  corrected  by  the  methods  that  follow. 

Wing-fuselage.-  The  corrections  for  wing-fuselage  contributions  are 

(1)  Correction  for  wing  -  From  figure  8  (based  on  reference  25) 
the  theoretical  values  of  C^p^Cp,  for  the  design  under  consideration 

and  for  the  design  for  which  data  are  available  can  be  determined.  The 
difference  between  these  two  theoretical  values  can  then  be  added  (with 
proper  regard  for  sign)  to  the  experimental  data.  Consideration  should 
be  given  to  scale  effect,  airfoil  section,  and  surface  roughness  on  the 
value  of  CZq  for  highly  swept  wings.  The  lift  coefficient  at  which 
the  experimental  variation  of  C^p  with  lift  coefficient  departs  from 

theory  is  greatest  at  high  Reynolds  numbers  and  for  smooth  wings  with 
round  leading  edges.  For  wings  with  rough  surfaces  or  sharp  leading 
edges  the  effects  of  Reynolds  number  on  Czp  are  usually  small  and  low- 

scale  wind  tunnel  data  can  be  used.  For  airplanes  having  very  smooth 
sweptback  wings  with  rounded  leading  edges,  however,  some  correction 
should  be  made  for  scale  effect  when  estimations  are  made  from  low-^ 
scale  wind-tunnel  data.  Since  no  rational  method  has  been  developed 
for  making  such  corrections  it  is  suggested  that,  for  lift  coefficients 
higher  than  that  at  which  the  experimental  data  departs  from  the  theory, 
an  average  of  the  theoretical  and  low-scale  experimental  values  be  used. 
Conservative  dynamic  stability  results  will  usually  be  obtained  if  the 
uncorrected  theoretical  values  of  are  used  because  these  values 

are  ordinarily  greater  (more  negative)  than  measured  values  and  because 
the  larger  negative  values  of  Czp  usually  tend  to  decrease  the  dynamic 
J.ateral  stability. 
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(2)  Correction  for  wing  dihedral  -  The  effect  of  dihedral  on 
is  treated  in  references  29^  39^  51^  58^  66^  and  79*  Correction  for 
the  difference  in  dihedral  between  the  two  designs  can  be  made  by 
multiplying  the  incremental  geometric  dihedral  angle  (in  degrees)  by 
the  factor  obtained  from  figure  9*  A  plot  of  against 

aspect  ratio  for  taper  ratios  of  1.0,  0.5  and  0.25  (obtained  from  refer¬ 
ences  58  and  66)  and  a  formula  from  reference  50  ^*0^  correcting  for 
sweep  are  presented  in  the  upper  portion  of  figure  9*  Ide  lower  chart 
and  formula  in  figure  9  (developed  from  reference  66)  should  be  used  in 
addition  to  the  upper  chart  and  formula  of  figure  9  estimate  the 
values  of  for  the  case  of  a  wing  with  partial-span  dihedral. 


Although  this  chart  and  formula  apply  directly  only  to  wings  with  one 
dihedral  break  they  can  be  used  to  estimate  the  C^p^  for  wings  with 

two  or  more  dihedral  breaks  by  the  method  described  in  reference  66. 
The  effect  of  drooped  wing  tips  and  of  wing-tip  end-plates  on 

should  be  determined  by  experimental  data  since  no  reliable  estimation 
procedure  for  these  effects  is  available. 


(3)  Correction  for  wing-fuselage  interference  -  Although  the  con¬ 
tribution  of  the  fuselage  alone  to  C^p  is  usually  negligible,  the 

interference  between  the  wing  and  fuselage  can  greatly  alter  the  value 
of  Czp  of  the  wing.  This  interference  is  such  that  a  high  location 

of  the  wing  on  the  fuselage  gives  more  positive  effective  dihedral 
(higher  -CZp)  and  a  low  wing  location  gives  less  positive  dihedral 
than  a  midwing  position.  This  effect  is  treated  theoretically  in 
reference  67  and  has  been  studied  experimentally  in  references  38  to  42, 
The  following  simplified  expression  for  estimating  the  increment  in  Czp 

caused  by  wing-fuselage  interference  has  been  developed  from  the  rela¬ 
tionships  presented  in  reference  67  and  in  other  sources: 


This  expression  has  been  found  to  give  reasonably  good  agreement  with 
experimental  data  for  a  variety  of  configurations.  It  is  suggested 
that  values  of  ACzp  be  calculated  from  this  equation  for  both  the 

design  under  consideration  and  for  the  design  for  which  force-test  data 
are  available.  The  difference  between  these  values  can  then  be  added 
(with  the  proper  regard  for  sign)  to  the  force -test  data. 

Vertical  tail.-  The  value  of  determined  from  force-test 

data  on  a  similar  design  can  be  corrected  as  follows  to  obtain  ^Zp^o^^p 
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for  the  design  under  consideration: 


Ptail 


design 


=  Cz 


Ptail 


\  Ptall  ^/design 

data 


3tail  ^/data 


The  results  of  reference  35  indicate  that  also  be  affected 

by  the  location  of  the  horizontal  tall  with  respect  to  the  vertical 
tail.  If  the  two  designs  have  approximately  the  same  horizontal  tail 
size  and  location,  however,  this  effect  can  be  neglected. 


The  value  of 

following  empirical  formula: 


for  a  V-tail  can  be  estimated  from  the 


^  _ 

r i)  si'n'^K-tail-^  W-tail  4 


'  \  _  /  \  JLb  sin  r  V^-tail"^  ^^V-tail  yjdesign 

^  ^V-tail /^0 3 \  ^V-tail)  (ig^ta  ^^PV-tail/  i  \ 

b  data 


(19) 

where  by-tail  is  the  developed  (not  projected)  span  of  the  V-tail, 
zv-tail  is  the  vertical  distance  from  the  center  of  gravity  to  the 
chord  of  the  V-tail  (positive  up,  and  F  is  the  dihedral  angle  of 
the.  V-tail.  More  information  on  V-tails  can  be  found  in  references  30, 
6l,  and  62. 

In  the  case  of  a  vertical  tail  located  on  the  wing,  there  is,  in 
addition  to  the  incremental  produced  by  the  tail  lateral  force, 

an  incremental  produced  by  the  interference  effect  of  the  vertical 

tail  on  the  wing.  Since  this  interference  effect  varies  greatly  with 
spanwise  and  vertical  position  of  the  tail,  it  should  be  determined 
from  force  tests.  Usually  the  in'cerference  is  such  that  a  vertical  tail 
above  the  wing  gives  a  negative  _ncrement  of  C^p  (positive  effective 

dihedral)  and  one  below  the  wing  gives  a  positive  increment  of  0^^. 

In  general,  the  largest  interference  effects  are  obtained  with  vertical 
tails  at  or  near  the  wing  tips. 
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THE  ROLLING  DERIVATIVES  Cn^,  Cj  ,  Cy 

Jr  ir 


The  wing  and  vertical  tail  are  the  only  airplane  components  that 
contribute  appreciably  to  the  yawing  moment  due  to  rolling  derivative 
Cn-n*  contributions  of  the  fuselage  and  horizontal  tail  can  usually 

sr 

be  neglected. 

Wing.  -  The  contribution  of  the  wing  to  Cnp  can  be  estimated  from 

the  formula  and  charts  of  figure  10  which  were  taken  from  reference  86. 
Although  these  charts  apply  strictly  only  to  wings  having  a  taper  ratio 
of  1.0,  experimental  data  have  indicated  that  they  will  also  provide 
fairly  good  estimates  for  taper  ratios  of  0.50,  0.25  and  0.  In  the 
estimation  formula 


the  value  of  should  be  determined,  if  possible,  from  force -test 

data  obtained  at  high  Reynolds  number  on  the  wing  under  consideration, 
since  low  Reynolds  number  data  might  indicate  values  of  that  are 

too  large.  For  the  case  of  smooth  wings  with  a  large  leading  edge 
radius  and  low  or  moderate  sweep,  it  is  suggested  that  for  the 

airplane  be  assumed  to  be  zero  at  all  lift  coefficients  up  to  the  stall. 
This  assuvaption  will  result  in  larger  negative  values  of  C^p  than 

would  be  estimated  from  low  Reynolds  number  data  on  and  con¬ 

sequently  should  lead  to  conservative  dynamic  stability  results  since 
an  increase  in  Cnp  in  the  negative  direction  has  been  found  to  cause 

a  reduction  in  dynamic  stability.  The  value  of  for  highly  swept 

wings  is  often  very  large  at  high  lift  coefficients,  especially  for 
wings  with  rough  surfaces,  sharp  leading  edges,  or  triangular  plan  form. 
For  these  cases,  values  of  determined  even  from  low  Reynolds 

number  data  might  lead  to  reasonably  good  estimates  of  Cup.  In  all 

these  cases,  however,  high-scale  drag  data  should  be  used  whenever  it 
is  available. 

Effect  of  high-lift  devices.-  The  principal  effect  of  leading-edge 
high-lift  devices  is  to  extend  to  a  higher  lift  coefficient  the  linear 
variation  of  Cn-rs  with  lift  coefficient.  The  formula  and  charts  of 

figure  10  are  directly  applicable  to  this  case.  The  effect  of 
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trailing-edge  high-lift  devices  is  not  so  straightforward,  hut  experi¬ 
mental  data  have  indicated  that  the  formula  and  charts  of  figure  10 
also  give  reasonably  good  estimates  in  this  case. 

Vertical  tail.-  The  contribution  of  an  isolated  vertical  tail 
surface  to  Cnp  can  be  estimated  by  the  following  approximate  formula 

which  has  also  been  commonly  used  to  estimate  of  a  complete 

airplane : 


=  -2  i  ^  CY. 


The  values  of  Cy 


Ptail 


'""Ptail  1  -"Ptail 

should  be  determined  from  force -test  data  as 


(21) 


*  . 

previously  discussed.  Instead  of  the  geometric  tail  length  ^/b,  it 
will  usually  be  better  to  use  the  effective  tail  length 

as  determined  by  force -test  data.  Formula  (21)  then  becomes 


'Xail  ^  H^f^ail 


(21a) 


In  the  case  of  the  conventionally  located  vertical  tall  surface,  how¬ 
ever,  the  rolling  wing  produces  a  sldewash  at  the  tall  which  greatly 


alters  the  tail  contribution  to 
of 


-n. 


This  sidewash  causes  the  values 


C„  to  be  much  more  negative  than  is  indicated  by  formula  (21). 

Ptail 


This  effect  is  discussed  more  fully  in  reference  36  in  which  is  also 
presented  a  method  for  estimating  the  sidewash.  Some  preliminary  theo¬ 
retical  studies  have  indicated  that  the  effect  of  the  sidewash  on 
Cn  varies  considerably  with  tail  size  and  tail  location  and  to 

^Ptail  ^ 

some  extent  with  wing  plan  form.  A  comprehensive  experimental  verifi¬ 
cation  of  this  theory  is  planned  but  as  yet  only  a  few  scattered  checks 
have  been  obtained.  For  the  case  of  the  conventionally  located  vertical 
tail  surface,  the  following  formula  has  been  found  to  give  estimates 
of  Cru  that  are  in  fairly  good  agreement  with  experimental  data: 

^tail 


or 


^tail 


^  i  i 

^  b  b 


^"Ptail 


(22) 


^^Ptail 


(22a) 
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This  formula  is  based  on  the  assumption  that  Cni 


is  zero  at  0°  angle 


^Ptail 

of  attack  and  varies  with  angle  of  attack  in  the  same  manner  as  indicated 
by  formula  (21).  Formula  (22)  or  the  method  of  reference  36  can  be  used 
satisfactorily  for  first  approximations  of  for  most  configura¬ 

tions  with  conventionally  located  vertical  tails.  For  more  accurate 
estimates,  especially  for  configurations  having  an  unusual  tail  size  or 
tail  location,  experimental  data  should  be  used. 


For  wings  of  triangular  plan  form  with  vertical  tails  either 
directly  above  or  above  and  slightly  behind  the  wing,  experimental  data 
have  indicated  that  neither  formula  (21)  nor  formula  (22)  gives  an 
accurate  estimate  of  an  average  of  the  values  obtained 

by  the  two  formulas  provides  a  fairly  good  estimate. 


It  is  obvious  that  these  methods  of  estimating  Cnp  are  only 

approximate  and  are  open  to  question  in  many  cases.  Experimental  and 
theoretical  studies  are  currently  being  made  to  provide  better  methods 
of  estimating  ^ap^g^^^p  and,  when  these  methods  become  available,  the 

approximate  methods  presented  herein  should  be  discarded.  At  the  present 
time,  however,  formula  (22)  and  reference  36  will  usually  provide  much 
more  accurate  estimates  of  ^np^^^^  than  formula  (21)  which  has  been 

in  common  use  up  until  this  time. 


Wing -fuselage . -  Most  of  the  rolling  moment  due  to  rolling  (damping- 
in-roll  derivative)  of  an  airplane  is  produced  by  the  wing.  The 

effect  of  the  fuselage  can  be  neglected  unless  the  ratio  of  the  diameter 
of  the  fuselage  to  the  wing  span  is  relatively  large  (greater  than 
about  0.3).  For  large  values  of  this  ratio,  the  value  of  will  be 

smaller  than  that  for  the  wing  alone  by  an  amount  that  can  be  estimated 
from  a  consideration  of  the  area  and  lateral  center  of  pressure  of  the 
wing  area  included  within  the  fuselage.  (See  references  103,  108, 
and  112.) 

Wing. -  The  damping  in  roll  of  wings  has  been  the  subject  of  many 
experimental  and  theoretical  investigations.  (See  references  on  CZp 
in  table  II.)  As  a  result,  some  methods  of  estimating  CZp  have  been 

developed  which  have  been  found  to  give  reasonably  good  agreement  with 
experimental  results.  The  method  presented  in  reference  79  appears  to 
give  sufficiently  accurate  estimates  of  Ci  for  zero  lift.  This 
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method  is  extended  in  reference  89  to  permit  the  estimation  of  C^p 

over  the  normal  flight  range  of  lift  coefficient.  Estimation  charts 
and  formulas  from  reference  89  are  presented  in  figure  11. 

Hlgh-lift  devices. -  Experimental  data  have  indicated  that  the 
damping  in  roll  of  wings  at  low  and  moderate  lift  coefficients  is  not 
greatly  affected  hy  the  addition  of  hi^-lift  devices  such  as  trailing- 
edge  flaps^  leading-edge  flaps,  slats,  and  slots.  The  principal  effect 
of  such  devices  is  to  increase  the  lift  coefficient  at  which  the  sharp 
decrease  in  CZp  occurs.  The  charts  and  formulas  of  figure  11  can  be 

used  to  estimate  the  Cl  of  wings  with  either  full-span  or  partial- 

xr 

span  nigh-lift  devices  with  fair  accuracy  despite  the  fact  that  the 
method  is  not  strictly  applicable  to  partial-span  high-lift  devices. 
(See  reference  89.) 

Wing-tip  fuel  tanks.-  The  use  of  wing -tip  fuel  tanks  usually 
increases  the  damping  in  roll  of  the  wing.  The  experimental  data  of 
reference  9I  unswept  wings  indicate  that  the  magnitude  of  the 
increase  varies  with  angle  of  attack  and  depends  upon  the  wing  taper 
ratio  and  on  the  size  and  location  of  the  tanks.  Unpublished  experi¬ 
mental  data  indicate  similar  effects  of  wing-tip  tanks  on  sweptback 
wings.  The  following  approximate  formula  for  estimating  the  increment 
in  CZp  produced  by  wing-tip  tanks  at  low  lift  coefficients  is  based 

on  the  limited  amount  of  available  experimental  data  and  should  not  be 
expected  to  yield  very  close  quantitative  estimates: 


tanks 


tanks  off 


Maximum  tank  diameter\  ^  ^ 


-Wing  span 


(23) 


where,  for  symmetrically  mounted  tip  tanks, 

Kt  =  6 

for  tanks  mounted  below  the  wing  tip  or  forward  on  the  wing  tip, 

Kt  =  3 

and  for  pylon -mounted  tip  tanks, 

Kj  =  1 

Experimental  data  for  both  unswept  and  swept  wings  indicate  that 

I  ACi  ]  .  usually  becomes  smaller  with  increasing  angle  of  attack 

\  y  x>  anx  s 

and,  in  some  cases,  actually  reverses  sign  at  high  angles  of  attack  so 
that  the  tanks  are  decreasing  rather  than  increasing  the  damping  in  roll. 
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The  data  of  reference  9^  can  be  used  to  obtain  an  approximate  estimate 
of  the  effect  of  angle  of  attack  for  unswept  wings. 

Tail  surfaces.-  The  contribution  uo  C^p  of  conventional  type 

horizontal  and  vertical  tail  surfaces  is  usually  very  small  and^  in 
most  cases^  negligible.  When  an  airplane  rolls^  the  wing  produces  a 
rotation  of  flow  at  the  tail  surfaces  which  reduces  the  already  small 
damping  moments  of  the  isolated  surfaces,  except  in  the  case  of  the 
vertical  tail  at  high  angles  of  attack  where  the  tail  center  of  pressure 
is  below  the  center  of  gravity. 

The  contribution  of  an  extremely  large  horizontal  tail  to  CZp 

might  not  be  negligible  and  can  be  estimated  by  multiplying  the  value 
of  C7  for  the  particular  tail  plan  form  obtained  from  the  charts  and 

in  which  the  fac¬ 
tor  0.5  is  included  to  account  for  the  rotation  of  flow  produced  by  the 
wing. 


formulas  of  figure  11  by  the  factor  0.5 


S  \b 


The  contribution  of  an  isolated  vertical  tail  surface  to 
given  by  the  following  approximate  formula: 


Czp  is 

(2^) 


As  in  the  case  of  formula  can  be  modified  to  provide  an 

approximate  correction  for  the  effect  of  the  wing  on  the  damping  in 
roll  of  conventionally  located  vertical  tail  surfaces: 


*^^Ptail  Kb) 

An  analysis  of  this  expression  i 

negligible  at  low  and  moderate  a 
but  that  it  might  be  fairly  important  at  very  high  angles  of  attack 
where  z/b  is  a  large  negative  value.  As  in  the  case  of  Cnp,  experi¬ 
mental  data  indicate  that,  for  a  vertical  tail  located  either  directly 
above  or  slightly  behind  a  wing  of  triangular  plan  form,  the  value  of 
can  be  estimated  with  better  accuracy  by  an  average  of  . 

i^tail 

formulas  (24)  and  (25)  than  by  formula  (25)  alone.  For  conventional 
tail  arrangements,  however,  formula  25  gives  better  correlation  with 
experimental  data. 


(t) 


0=0 


tall 


(25) 


idicates  that  the  value  of  C?  is 

^Ptail 

igles  of  attack  where  z/b  is  positive 
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Cy 

Wing.-  The  following  formula  for  the  derivative  Cyp  (lateral 

force  due  to  rolling)  from  reference  86  is  based  on  experimental  data 
and  is  the  same  as  that  presented  in  reference  25  except  for  an  addi¬ 
tional  correction  to  account  for  tip  suction: 


Cl 


A  +  cos  A 
A  +  4  cos^^ 


tan  A  + 


(26) 


The  data  of  reference  86  show  that  this  formula  applies  only  for  lift 

Cl^ 

coefficients  below  that  at  which  the  drag  factor  Cj)  -  begins  to 

Increase.  At  hi^er  lift  coefficients  the  experimental  data  Indicate 
smaller  values  of  Cyp  than  given  by  formula  (26).  For  these  cases  an 

approximation  of  the  value  of  Cyp  can  be  obtained  from  the  experi¬ 
mental  data  of  reference  86.  As  in  the  case  of  Cnp,  the  break  in  the 
variation  of  Cyp  with  lift  coefficient  should  be  expected  to  occur 

at  lower  lift  coefficients  for  wings  having  sharp  leading  edges  or 
rough  surfaces  and  for  wings  tested  at  low  Reynolds  numbers. 


Vertical  tail.-  The  discussion  concerning 
is  also  applicable  to  The  value  of 

tail  surface  is  given  by  the  formula: 


and  C7r, 

Ptail 

for  an  Isolated 


CYn 


tail 


(l)cYp 


tail 


(27) 


This  formula  can  be  modified  as  follows  to  account  approximately  for  the 
effects  of  wing  sidewash  in  the  case  of  a  conventionally  located  vertical 
tail: 


tail 


(28) 


An  average  of  formulas  (27)  and  (28)  can  be  used  for  tails  located 
either  directly  above  or  above  and  slightly  behind  the  wing. 
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THE  YAWING  DERIVATIVES  AND  Cy^ 


Cnr 


Wing-fuselage.-  In  the  past,  the  contribution  of  the  wing -fuselage 
combination  to  yawing  moment  due  to  yawing  (damping  in  yaw)  derivative 
Crij.  has  usually  been  found  to  be  small  compared  to  the  contribution  of 

the  vertical  tail.  The  fuselage  contribution  to  the  damping  in  yaw 
depends,  of  course,  on  the  relative  size  of  the  fuselage  and  wing.  In 
the  past,  the  relative  size  of  these  components  has  generally  been  such 
that  the  fuselage  contribution  could  be  neglected.  (See  references  82 
and  83.)  For  some  recent  designs  which  have  a  large  i’uselage  relative 
to  the  wing,  however,  the  fuselage  contribution  to  Cnp  is  important. 

In  the  case  of  fuselages  having  flat  sides  or  having  a  flattened  cross 
section  with  the  major  axis  vertical  the  fuselage  contribution  may  also 
be  important  and  some  fuselage  contribution  to  Cn.j.  should  be  assumed, 

especially  at  high  angles  of  attack.  On  the  other  hand,  experimental 
data  have  shown  that  a  flattened  cross-section  fuselage  with  the  major 
axis  horizontal  can  have  negative  damping  in  yaw  at  moderate  and  high 
angles  of  attack. 

The  contribution  of  the  wing  to  Cnr  can  be  estimated  from  the 
formula  and  charts  of  figure  12  which  were  taken  from  reference  25. 
Values  of  for  the  wing  should  be  estimated  from  force-test  data. 

For  values  of  x/c  greatly  different  from  zero^  the  charts  of  refer¬ 
ence  25  can  he  used.  The  formula  and  charts  of  figure  12  are  not  con¬ 
sidered  reliable  at  high  angles  of  attack,  especially  for  swept  wings. 
The  use  of  experimental  data  from  the  references  on  Cn^.  listed  in 

table  II  is  recommended  in  this  case. 

The  effect  of  partial-span  inboard  flaps  on  Cn^.  can  usually  be 

neglected.  (See  reference  82.)  The  effect  of  full-span  trailing-edge 
or  leading-edge  hi^-lift  devices  can  be  estimated  satisfactorily  from 
the  formula  and  charts  of  figure  12. '  Values  of  CDq  in  this  case  are, 

of  course,  for  the  wing  with  the  high-lift  device  installed. 

Vertical  tall.-  The  contribution  of  a  conventional -type  vertical 
tail  to  can  be  estimated  from  the  formula 


^""^tall  ^(h)  ^^Ptail 


(29) 
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or,  with  the  effective  tail  length  “^'^3-ta,ii/*^^Ptail 
the  geometric  tail  length  z/h. 


'^Ptall. 

'^Ptail 


(29a) 


The  experimental  values  for  Cm-  presented  in  reference  82  for 

■^tail 

power-on  or  flap-down  configurations  are  30  to  40  percent  greater  than 

values  predicted  by  formulas  (29)  or  (29a).  These  differences  are 

attributed  to  lag  of  sidewash  effects  in  the  free-oscillation  tests 

used  in  measuring  Cm-.  In  estimations  of  C^^  for  stability  cal- 

^  tail 

culations^  similar  lag  of  sidewash  effects  should  be  assumed  if  the 
oscillatory  mode  is  of  primary  importance  but  no  lag  of  sidewash  should 
be  assumed  if  the  aperiodic  mode  is  most  important. 

Methods  for  estimating  the  wing-tip  vertical  tails 

are  presented  in  references  70  and  82. 


The  wing  and  vertical  tail  are  the  only  airplane  components  that 
contribute  appreciably  to  rolling-moment -due -to -yawing  derivative  C^^ 

of  an  airplane.  The  contributions  of  the  fuselage  and  horizontal  tail 
can  usually  be  neglected.  A  semiempirical  method  for  estimating  Ci^ 

is  presented  in  reference  85.  This  method  involves  the  use  of  experi¬ 
mental  data  on  the  parameter  C^p  to  correct  the  theoretical  values  of 

C^j^wing  reference  25  and  to  estimate  the  value  of 

Wing.  -  The  formula  of  reference  85  and  the  charts  of  Ci^jCij  from 


reference  25  for  estimating  C^ 


are  given  in  figure  I3 


■^wing 

values  of  to  be  used  in  the  charts  can  be  obtained  from  fig¬ 

ure  8.  For  taper  ratios  less  than  0.25^  values  of  C^^^/Cl  and 
for  a  taper  ratio  of  0.25  can  be  used.  The  value  of 


the  formula  should  be  the  same  as  the  value  of  C^ 


n  UDCU.  J-ii. 

Pexp 

estimated  from 


wing 


experimental  data  by  the  method  indicated  in  the  section  on  In 

the  case  of  however,  (unlike  the  case  of  Czp)  conservative 
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dynamic  stability  results  will  usually  be  obtained  if  the  smaller  values 
of  the  derivative  (based  on  low-scale  experimental  data)  are  used 
instead  of  the  larger  (theoretical)  values.  This  difference  is  a  result 
of  the  fact  that  either  an  increase  in  the  normally  negative  value  of 
CZp  or  a  decrease  in  the  normally  positive  value  of  can  cause 

reduction  in  dynamic  stability.  As  pointed  out  in  reference  85  the 
estimation  procedure  shown  in  figure  I3  appears  to  account  satisfactorily 
for  the  effects  of  high-lift  devices,  wing,  dihedral,  and  airfoil  section. 

Vertical  tail.-  The  contribution  of  the  vertical  tall  to  Ci^  is 
usually  estimated  by  the  formula 

‘^^I'tail  "  ■^(F)(s)^Y3tail  ^30) 

where  is  preferably  obtained  from  force-test  data.  When 

experimental  data  on  available,  the  following  fonnula  from 

reference  85  can  be  used  and  will  probably  be  more  reliable  than  equa¬ 
tion  (30)  because  it  takes  into  account  any  interference  effects  that 
mi^t  cause  the  effective  vertical  location  of  the  center  of  pressure  of 
the  tail  to  be  different  from  the  location  determined  by  geometrical 
procedures: 


Cl 


^tail 


^Ptail 


(31) 


or  with  the  effective  tail  length  ”*^'^Ptall/'^^3tail  substituted  for  the 
geometric  tall  length  Z/b, 


Cl  =  2I 

^tail 


Ptail 


.CYp 


Cl, 


tail^ 


’tail 


(31a) 


Cy 


r 


Wing. -  The  theory  of  reference  25  gives  values  of  the  derivative 
CYp  (lateral  force  due  to  yaving)  for  the  wing  for  a  taper  ratio  of  1.0. 

The  experimental  data  of  references  25  and  59  indicate  that  this  theory 
is  inadequate  for  making  reliable  estimates  of  Cy^  .  It  is  recom- 

^wing 

mended  therefore  that  the  experimental  data  given  in  references  25,  58 
59 >  f‘'ind  60  be  used  in  making  estimates  of  ^Yr^in.  *  ^  ^ 


The  effects  of  Mach  number  on  the  lateral  stability  derivatives 
have  been  treated  theoretically  in  many  investigations  (see  table  II ) 
but  very  little  experimental  data  have  been  obtained  to  verify  this 
theoretical  work.  Moreover,  only  a  small  part  of  this  experimental 
work  has  been  covered  in  published  reports  (reference  111)  because  most 
of  it  is  classified  at  the  present  time.  It  appears,  therefore,  that 
estimates  of  the  lateral -stability  derivatives  for  the  time  being  will 
have  to  be  based  largely  on  theoretical  work. 

The  effects  of  Mach  number  on  the  stability  derivatives  can  be 
usually  considered  negligible  for  all  airplane  components  except  the 
wing  and  vertical  tail.  For  the  low-lift-coefficient  condition  in  the 
case  of  many  high-speed  airplanes,  the  vertical  tail  contributes  more 
than  the  wing  to  all  the  stability  derivatives  except  CZp-  For  this 

reason,  in  calculations  for  transonic  or  supersonic  speed  conditions  it 
is  especially  important  to  know  the  effects  of  Mach  number  on  the 
vertical -tail  lift-curve  slope  or  • 

Wing.-  The  effects  of  compressibility  on  the  subsonic  stability 
derivatives  of  the  wing  can  be  estimated  by  the  formulas  of  reference  26. 
The  values  of  the  supersonic  stability  derivatives  for  some  wing  plan 
forms  can  be  estimated  by  the  references  tabulated  in  table  II.  In  this 
table  the  derivatives  are  grouped  according  to  the  type  of  wing  plan 
form  and  to  the  particular  derivatives  covered.  A  helpful  summary  and 
discussion  of  the  effects  of  Mach  number  on  the  derivatives  for  several 
different  wing  plan  forms  is  presented  in  reference  103.  A  summary  of 
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the  theoretical  lift -curve  slope,  dajnping  in  roll,  and  center-of- 
pressure  characteristics  of  various  wing  plan  forms  is  presented  in 
reference  107 .  Iti  the  cases  in  which  the  theory  shows  large  or  abrupt 
changes  in  a  stability  derivative  with  changes  in  Mach  number  (for 
example,  fig.  10  of  reference  103)  special  care  should  be  taken  in 
estimating  the  derivative  in  that  particular  Mach  number  range.  The 
abrupt  changes  should  be  smoothed  or  faired  out  in  a  manner  similar  to 
that  suggested  in  the  following  section  for  estimating 

In  some  cases,  experimental  data  for  supersonic  speeds  will  be 
available  on  the  sideslip  derivatives  and  on  the  damping -in-roll  deriva¬ 
tive  Cip.  In  such  cases  the  experimental  data  should  be  used  in  pref¬ 
erence  to  the  theory.  Some  experimental  results  have  indicated  that  the 
effect  of  the  vertical  location  of  the  wing  on  the  fuselage  on  the 
derivative  might  be  greatly  different  at  supersonic  speeds  from 

that  at  subsonic  speeds.  Since  no  methods  are  presently  available  for 
estimating  this  effect  for  the  supersonic  case,  it  appears  that,  at 
least  in  the  case  of  high -wing  and  low-wing  designs,  force -test  data 
are  necessary  for  obtaining  an  accurate  estimate  of  C^p. 


Vertical  tail.-  The  sideslip  derivatives  produced  by  the  vertical 
tail  at  transonic  and  supersonic  speeds  can  be  estimated  theoretically 
but  should  be  obtained  from  force-test  data  whenever  possible.  These 
sideslip  derivatives  can  be  used  to  estimate  the  tail  contributions  to 
the  other  derivatives  as  pointed  out  previously.  In  estimates  of  the 

value  of  Cyq  for  transonic  and  supersonic  speeds,  corrections  must 

*^tail 

be  made  for  the  effect  of  Mach  number  on  the  lift'-curve  slope  of  the 
tail,  and  these  corrections  should  account  for  any  differences  in  the 
end-plate  effect  of  the  horizontal  tail  on  the  vertical  tail. 


For  Mach  numbers  below  about  0.8  or  0.9  and  above  about  1.6  or  1.8 
the  effect  of  Mach  number  on  the  lift-curve  slope  of  the  vertical  tail 
can  be  estimated  satisfactorily  from  the  theoretical  values  of  refer¬ 
ences  26,  3^^  and  IO7.  Since  experimental  data  indicate  that  theoreti¬ 
cal  values  of  lift-curve  slope  are  usually  too  high  for  Mach  numbers 
from  about  0.8  or  0.9  to  about  1.6  or  1.8,  the  empirically  determined 
fairings  shown  in  figure  ik  are  recommended  for  use  as  a  guide  in  the 
use  of  the  theory  to  obtain  approximate  estimates  in  this  Mach  number 
range  when  force -test  data  are  not  avilable. 


Experimental  data  have  indicated  that  for  vertical -tail  configura¬ 
tions  which  have  a  tail  length  (distance  from  the  center  of  gravity  to 
the  tail  center  of  pressure)  that  is  relatively  short  in  terms  of  tail 
chords,  the  rearward  shift  of  the  tail  center  of  pressure  at  supersonic 
speeds  can  cause  an  appreciable  increase  in  the  tail  length  and 


k& 
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consequently  an  appreciable  Increase  in  the  magnitude  of  some  of  the 
tail  derivatives.  Theoretical  center-of -pressure  positions  for  various 
plan  forms  at  supersonic  speeds  are  given  in  reference  107 . 


EFFECTS  OF  POWER 


On  the  basis  of  existing  Information,  the  effects  of  power  on  the 
lateral  stability  derivatives  appear  to  be  negligible  in  the  case  of 
jet-propelled  airplanes  but  these  effects  are  often  very  large  in  the 
case  of  single-engine  propeller -driven  airplanes.  Methods  are  available 
for  estimating  some  of  these  power  effects  but  in  most  cases  experimental 
data  are  necessary  for  making  a  satisfactory  estimate.  The  effects  of 
pow6r  can  be  broken  dovn  inbo  t-wo  general  classes; 

(1)  The  effects  of  the  lateral  force  produced  by  the  propeller 

itself 

(2)  The  effects  of  the  propeller  slipstream  on  the  wing,  fuselage, 

and  vertical  tail  of  the  airplane 

Effects  of  propeller  lateral  force.-  A  method  of  estimating  the 
propeller-lateral-force  derivative  Cyp  is  presented  in  reference  31 

which  is  based  on  the  work  of  references  32  and  33*  The  contribution 
of  the  propeller  lateral  force  to  the  other  stability  derivatives ^ can 
be  estimated  from  this  derivative  by  assuming  that  the  propeller  is 
effectively  a  vertical  tail  surface  and  by  using  the  expressions  for 
the  tail  contribution  to  the  various  derivatives  presented  in  the 
preceding  sections.  Some  experimental  data  on  the  effect  of  windmilling 
propeller  on  all  of  the  derivatives  are  presented  in  reference  65. 

Effects  of  propeller  slipstream.-  The  effects  of  propeller  slip¬ 
stream  on  the  lateral' stability  derivatives  are  usually  much ^ greater 
than  the  effects  of  propeller  lateral  force  in  the  case  of  single -engine 
tractor  airplanes.  The  slipstream  effects  on  the  wing,  the  fuselage, 
and  the  vertical  tail  can  be  considered  as  three  independent  effects. 


The  slipstream  effects  on  the  wing  can  usually  be  neglected  except 
for  the  derivatives  Czp  and  Ci^,  Experimental  data  showing  the 

decrease  in  effective  dihedral  with  power  for  single-engine  air¬ 

planes  are  presented  in  references  5^^  55^  5^,  and  80.  It  appears 
highly  desirable  to  determine  this  effect  of  power  experimentally 
because  interference  effects  make  accurate  estimations  of  the  effect 
very  difficult.  The  effect  of  the  slipstream  on  the  value  of 

cannot  be  estimated  from  the  data  on  described  in  the 
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section  on  In  fact^  this  procedure  would  probably  give  the  wrong 

sign  for  the  increment  of  ^^rving  “the  slipstream.  An 

approximation  of  this  increment  might  be  obtained  by  estimating  the  slip¬ 
stream  velocity  and  the  lateral  displacement  of  the  slipstream  caused 

by  yawing.  Usually  the  power  effects  on  C^o  and  Ci^  will  be 

^wing  ^ wing 

greatest  for  the  flap -extended  configuration. 

In  the  case  of  the  single -engine  airplane  the  effect  of  the  slip¬ 
stream  on  the  fuselage  is  usually  to  increase  negatively  the  values 
of  Cnp  and  Cyp.  (See  references  5^^  55^  5^^  71^  7^^  sjid  76.)  Since 

no  accurate  methods  of  estimating  these  slipstream  effects  on  Cnp 
and  Cyp  are  available^  it  is  necessary  to  determine  them  from  force- 
test  data. 

The  effects  of  the  slipstream  on  the  vertical  tail  are  often  very 
important  and  should  also  be  determined  from  experimental  data^  if 
possible.  The  increase  in  dynamic  pressure  at  the  tail  caused  by  the 
slipstream  is  treated  theoretically  in  reference  ll6  and  is  illustrated 
by  the  experimental  data  of  references  50^  5^^  55^  5^,  71^  7^^  76. 

The  experimental  data  of  reference  76  also  show  that  the  propeller  slip¬ 
stream  can  cause  a  destabilizing  sidewasji  at  the  tail  which  will  tend 
to  reduce  the  stabilizing  effect  of  the  increased  dynamic  pressure  at 
the  tail.  Since  these  data  indicate  that  slipstream  effects  on  the 
vertical  tail  vary  greatly  with  airplane  configuration  and  propeller 
arrangement  (single  or  dual  rotation)^  use  of  experimental  data  appears 
to  be  the  only  satisfactory  estimation  procedure  at  present. 

Suggested  estimation  procedure  for  power  effects. -  The  following 
procedure  is  suggested  for  estimating  power  effects.  Obtain  force- 
test  data  for  tail  off  and  tail  on.  Use  tail-on  data  directly  for  Cyp^ 

Cnp^  and  C^p.  Estimate  rolling  and  yawing  derivatives  as  follows: 

(1)  Estimate  Cyn  _  from  reference  3I  and  use  this  deriva- 

■^Ppropeller 

tive  and  proper  linear  dimensions  to  estimate  other  propeller  deriva¬ 
tives  (rolling  and  yawing  derivatives)  in  the  same  manner  as  tail 
derivatives . 

(2)  Subtract  tail-on  data  from  tail-off  data  to  get  values  of 

Cyptaii^  Cnptaii^  ^^Ptail  power-on  condition  and  use 

these  values  to  estimate  the  tail  contribution  to  the  other  derivatives. 

(3)  For  tail-off  values  of  rolling  and  yawing  derivatives,  use 
same  values  as  for  power-off  for  all  derivatives  except  C2j;..  Estimate 

C2p  as  suggested  in  preceding  section. 
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(A)  Add  the  values  obtained  in  steps  1,  2^  and  3  to  get  the 
rolling  and  yawing  derivatives  for  the  complete  airplane. 

INADEQUACIES  IN  PRESENT  INFORMATION  AND  METHODS 


In  the  course  of  summarizing  the  estimation  methods  for  the  various 
stability  derivatives,  the  need  for  much  additional  information  on  all 
the  derivatives  became  apparent.  In  particular,  information  is  needed 
to  aid  in  the  estimation  of  the  derivatives  in  the  transonic  and  super¬ 
sonic  speed  ranges.  Additional  work  also  needs  to  be  done  in  correlating 
and  analyzing  existing  subsonic  data  and  in  obtaining  new  experimental 
data  for  the  development  of  semiempirical  methods  of  estimating  the  sub¬ 
sonic  derivatives  without  resort  to  force -test  data.  Another  important 
need  is  for  full-scale  experimental  results  at  all  speeds  for  checking 
both  low-scale  data  and  the  existing  methods  of  estimating  derivatives. 
Details  of  the  need  for  additional  work  along  these  lines  are  discussed 
in  the  following  sections.  Studies  should  also  be  made  to  determine 
the  conditions  for  which  the  use  of  steady-state  stability  derivatives 
in  conventional  stability  equations  is  inadequate  and  to  determine 
satisfactory  methods  of  treating  such  conditions. 


Transonic  and  Supersonic  Speeds 

Additional  theoretical  work  is  needed  on  the  estimation  of  sta¬ 
bility  derivatives  in  the  transonic  and  supersonic  speed  ranges  to 
cover  the  range  of  wing  plan  forms  for  all  the  derivatives.  In  particu¬ 
lar,  more  work  is  needed  on  plan  forms  currently  under  consideration, 
such  as  wings  having  moderate  sweepback  and  taper.  This  need  is  illus¬ 
trated  by  table  II  which  indicates  that  very  little  material  is  available 
on  the  stability  derivatives  for  such  plan  forms  except,  perhaps,  for 
the  derivative  CZp.  It  appears  from  the  table  that  this  derivative 

and  the  triangular  plan  form  have,  in  the  past,  received  a  dispropor¬ 
tionate  share  of  attention,  probably  because  of  the  greater  ease  with 
which  they  could  be  treated  theoretically. 

The  greatest  need  for  work  on  stability  derivatives  at  the  present 
time  is  probably  in  the  measurement  of  the  derivatives  at  transonic  and 
supersonic  speeds.  Experimental  data  on  wings  are  urgently  needed  for 
checking  the  theoretical  work  and  for  use  in  the  development  of  empirical 
corrections  to  the  theory  wherever  necessary.  Such  corrections  are 
particularly  needed  for  fairing  out  abrupt  variations  of  the  derivatives 
with  Mach  number  and  for  fairing  through  the  Mach  number  range  for  which 
theory  predicts  infinite  values.  Examples  of  such  discontinuities  as 
indicated  by  theory  are  shown  in  figures  8  to  13  of  reference  103- 
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Since  experimental  data  obtained  at  supersonic  speeds  on  wing-fuselage 
combinations  and  on  complete  models  have  revealed  interference  effects 
that  are  different  from  those  obtained  at  subsonic  speeds^  it  appears 
highly  desirable  to  obtain  at  least  a  limited  amount  of  experimental 
data  at  transonic  and  supersonic  speeds  to  evaluate  these  interference 
effects.  For  example^  investigations  should  be  undertaken  to  determine 
the  effect  of  wing-fuselage  interference  on  the  derivative  Ci^  and 

the  end-plate  effect  of  the  horizontal  tail  on  the  lift-curve  slope  of 
the  vertical  tail. 

Most  of  the  experimental  data  on  stability  derivatives  at  transonic 
and  supersonic  speeds  will  of  necessity  be  obtained  at  Reynolds  numbers 
considerably  less  than  full-scale  values  and  under  test  conditions  which 
might  render  the  results  open  to  question  in  some  cases.  Full-scale 
checks  in  flight  of  the  low- scale  data  and  of  the  estimation  methods 
therefore  appear  to  be  desirable.  Consequently  the  methods  of  measuring 
stability  derivatives  in  flight  now  being  developed  by  the  Cornell 
Aeronautical  Laboratory^  the  Massachusetts  Institute  of  Technology^  and 
the  NACA  should  be  extended  to  transonic  and  supersonic  speeds  when  the 
methods  appear  to  be  developed  to  a  satisfactory  degree  of  reliability 
for  the  subsonic  case.  Some  preliminary  considerations  involved  in  the 
use  of  these  flight  techniques  are  discussed  in  references  II7  to  120. 


Subsonic  Speeds 

The  methods  presented  in  this  paper  for  estimating  the  stability 
derivatives  at  subsonic  speeds  depend  either  directly  or  indirectly  on 
the  use  of  force -test  data.  These  methods  are  probably  more  reliable 
than  methods  which  do  not  involve  the  use  of  force -test  data  on  the 
particular  design  under  consideration  or  on  a  similar  design.  Methods 
which  do  not  rely  on  such  data  are  desirable  in  some  cases ^  however^ 
because  the  necessary  data  will  not  always  be  available. 

In  the  case  of  sideslip  derivatives,  empirical  methods  can  probably 
be  developed  largely  from  existing  information.  In  some  cases  it  will 
be  necessary  to  augment  the  existing  information  with  new  results  since 
much  of  the  available  force -test  data  were  not  obtained  in  a  manner  that 
would  make  the  data  readily  usable  for  developing  general  estimation 
procedures . 

In  the  case  of  rolling  and  yawing  derivatives,  considerably  less 
information  is  available  than  in  the  case  of  the  sideslip  derivatives. 
Most  of  the  infonnation  now  available  was  obtained  in  the  Langley 
stability  tunnel,  principally  on  wing  configurations  and  to  a  limited 
extent  on  complete  airplane  models  and  airplane  components  other  than 
the  wing.  Considerably  more  work  is  required,  especially  for  components 
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in  combination,  before  satisfactory  methods  can  be  developed  for 
estimating  rolling  and  yawing  derivatives  without  the  use  of  force-test 
data  on  the  particular  design  under  consideration  or  on  a  similar  design. 

In  discussing  the  work  necessary  for  developing  new  procedures  for 
estimating  the  stability  derivatives  without  the  use  of  force-test  data 
on  the  design  under  consideration  or  on  a  similar  design,  it  is  useful 
to  break  the  problem  down  into  two  parts:  (1)  effect  of  individual 
components  and  (2)  the  effect  of  Interference  of  the  components  on  each 
other. 

The  principctl  components  to  be  considered  are  the  fuselage^  wing^ 
vertical  tail^  and  propeller.  For  the  isolated  fuselage^  the  main 
problem  is  the  development  of  methods  for  the  estimation  of  and 

then^  perhaps^  of  Cn.^  sJ^d  CJyp-  For  the  isolated  wing^  the  main 

problem  is  to  estimate  the  derivatives  at  lift  coefficients  above  that 
at  which  separation  begins.  Such  estimations  can  be  made  with  reasonable 
accuracy  for  some  of  the  derivatives  by  existing  methods  which  make  use 
of  force-test  data^  but  the  development  of  methods  which  do  not  involve 
the  use  of  force-test  data  will  probably  be  very  difficult.  For  the 
isolated  vertical  tail^  the  problem  is  to  establish  the  effective  tail 
area  and  aspect  ratio  from  the  geometry  of  the  tail  so  that  the  lift*" 
curve  slope  ^or  oY  tail  can  be  calculated.  Solutions  to  this 

seemingly  simple  problem  have  in  the  past  become  involved  with  interfer¬ 
ence  effects  so  that^  as  yet^  no  reliable  methods  have  been  published 
for  estimating  Gyp  of  the  vertical  tail  from  its  geometry.  For  the 

isolated  propellers^  the  work  that  is  needed  at  present  is  a  systematic 
check  of  existing  methods  of  estimating  the  lateral  force  on  the 
propeller  to  determine  the  accuracy  of  these  methods. 

The  principal  interference  effects  to  be  considered  are  mutual 
interference  of  the  wing  and  fuselage ^  wing-fuselage  interference  on 
the  vertical  tail^  horizontal -tail  interference  on  the  vertical  tail^ 
propeller-slipstream  interference  on  the  wing^  fuselage^  and  vertical 
tail.  The  mutual -interference  effects  of  the  wing  and  fuselage  are 
probably  important  only  for  the  derivatives  Cnp^  and  A 

large  amount  of  experimental  data  is  available  for  the  sideslip  deriva- 
Uves  but  no  procedures  for  estimating  the  interference  effects  on  these 
derivatives  have  been  reported.  Wing-fuselage  interference  has  very 
important  effects  on  Gyp  of  the  vertical  tail,  and  consequently  on  all 

of  the  stability  derivatives  for  some  flight  conditions.  These  effects 
result  from  the  sidewash  and  change  in  dynamic  pressure  at  the  tail 
which  may  result  from  sideslipping,  rolling,  or  yawing.  Although  con¬ 
siderable  data  which  show  these  interference  effects  are  available, 
particularly  for  the  case  of  sideslipping,  no  reliable  methods  exist 
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for  estimating  the  interference  effects.  Horizontal -tail  interference 
also  has  an  important  effect  on  Cyp  of  the  vertical  tail  for  some 

horizontal-tail  positions.  Some  work  on  a  limited  number  of  configura¬ 
tions  has  been  done  toward  developing  methods  of  estimating  this  effect 
but  data  are  required  on  more  configurations  before  the  generally  appli¬ 
cable  methods  can  be  evolved.  The  propeller  slipstream  can  cause  impor¬ 
tant  effects  on  and  of  the  wing,  on  Cng  and  Cyg  of  the 

fuselage,  and  on  Cyg  of  the  tail  (and  consequently  on  the  tail  contri¬ 
bution  to  all  the  derivatives).  Some  data  are  available  for  the  effect 
of  the  slipstream  on  the  sideslip  derivatives  but,  because  of  the  com¬ 
plexity  of  this  problem,  considerable  additional  data  may  be  required 
before  a  satisfactory  method  of  estimating  the  slipstream  effects  can 
be  developed. 

As  mentioned  in  the  preceding  section,  full-scale  checks  of  low- 
scale  data  and  of  the  estimation  methods  are  desirable.  For  the  sub¬ 
sonic  case  some  of  the  checks  can  be  obtained  from  large-scale  wind- 
tunnel  tests  but  some  checks  in  full-scale  flight  tests  should  also  be 
obt'  -ined  when  the  various  methods  of  measuring  stability  derivatives  in 
flit,nt  have  been  developed  to  a  satisfactory  degree  of  accuracy. 


Langley  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Langley  Field,  Va.,  December  13,  1950 
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APPENDIX  A 


EQUATIONS  OF  MOTION 


The  dimensional  equations  for  the  lateral  motions  of  an  airplane 


mV  2  ^4  ^  SL  dt  SL  t  -  n 


SN  dci  ,  ,2  d^t  Sn  dt  SN 

'XZ  ^  ^  +  “kz  -  BF  dt  -  ^ 


V  -  N„  =  0 


H  It  "  ^  It  ■  If  it  ~  (Lift)(tan  7)t 


^  _  ax 

dt  ”  ^ 


If  equations  (Al)  and  (A2)  are  divided  by  — pV^Sb  and  equation  (A3)  is 
Ip  ^ 

divided  by  ^  equations  of  motion  may  be  expressed  in  the  con¬ 

ventional  nondimensional  form  in  which  they  have  generally  been  presented 
in  NACA  reports  (for  example^  see  reference  2): 


2,Kx"  ^  -  \  .  2.Kxz  g  -  I  C 


d^  d\lr  ^ 

^  “  2  IF  ■  ° 


pjifC  ^  i  n  +  2uKr7^  ^  ^  _  i  p  -  Q  3  -  C  =  0 

2^"%Z  .^2  2  ^np  ds  2  ds  ^nc  ^ 


1  CYr,  -  Cl0  +  2p 

2  ds  ^  ds 


CY33  -  Cy^  =  0 


I  CYr  g  -  CL(taii  7)^  +  g 
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In  order  to  convert  these  equations  into  a  form  which  will  reduce  the 
number  of  arithmetical  and  algebraic  steps  in  performing  stability  cal¬ 
culations^  equations  (a4)  are  multiplied  by  m/pSb  and  written  in  the 
following  form: 
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APPENDIX  B 


APPLICATION  OF  THE  LAPLACE  TRANSFORMATION  TO  CALCULATING  MOTIONS 


The  application  of  the  Laplace  transformation  to  the  calculation 
of  the  lateral  motions  of  airplanes  is  presented  in  order  to  illustrate 
the  development  of  the  equations  of  motion  in  the  form  in  which  they 
are  presented  in  the  present  paper.  This  work  is  similar  to  that  pre¬ 
sented  in  references  5  6.  In  fact,  it  follows  the  presentation  in 

reference  5  very  closely.  Reference  6  presents  a  brief  explanation  of 
the  Laplace  transformation  and  its  application  to  solution  of  the  equa¬ 
tions  of  motion  of  an  airplane.  This  paper  also  makes  reference  to 
detailed  explanations  of  the  Laplace  transformation.  In  cases  where 
modification  of  the  equations  presented  in  the  present  paper  are  neces¬ 
sary,  reference  should  be  made  to  these  texts  for  an  understanding  of 
the  mathematics  Involved.  Applying  the  Laplace  transforms 

L(l)  =  i  L(]D0)  =  -  00 

L(0)  =  0^  L(d20)  =  X%  -  X0O  -  (C0)o 


and  multiplying  each  of  the  equations  by  X  transforms  equations  (A5) 
from  appendix  A  to 


where 

ri  =  {X^  -  ^p^)!^0  ^c 

r2  =  (K2X^  -  npX)0o  +  {X^  -  nrX)io  +  K2A,(D0)o  +  X.(D’|f)o  +  Hc 
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Solving  equations  (Bl)  By  determinants  gives 


-ZpX 

^1 

KiX3  -  ZrX^ 

-npX 

x3  -  nj.X^ 

^3 

X^  -  yi.X^  -  ^(tan  y)X 

-ZpX, 

x3  -  ZpX^ 

KiX3  -  ZyX^ 

-n^A. 

K„x3  -  n  X^ 

x3  -  n  X^ 

(B 

2  p 

r 

X2  -  ypX 

■V  2  ^ 

x^  -  yj.^^  -  -^(tan  7)X 

which  may  be  expressed  as 

+  aj^X^  +  8^x3  +  a^X^  +  aj^X  +  a^ 

rX  ~  27  If  ^  ^  \ 

X  [AX  +  BX^  +  CX  +  DX  4-  Ej 

Similarly^  the  expressions  for  and  |3^  are 


ho^^  +  BjX^  +  b2X-^  +  1^3^^  +  ^4^  ^5 


AX^  +  BX^  +  CX^  +  DX  +  E 


(B3) 


CqX^  +  C]^x3  +  cgX^  +  CgX  +  cj^ 
X  +  Bx3  +  CX^  +  ,DX  +  E  ^ 


(b4) 


where  the  expressions  for  the  coefficients  in  equations  (B2)  to  (b4)  are 
given  in  terms  of  the  mass  and  aerodynamic  stability  derivatives  by 
equations  (l)  to  (4)  in  the  main  body  of  this  paper. 


In  order  to  obtain  the  actual  variables  from  the  transformed 
variables^  an  inverse  Laplace  transformation  must  be  applied.  The 
expressions  for  0^,  and  (3^  are  of  the  form  ^x/^X  ^x 
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and  v-^  are  polynomials,  the  degree  of  being  higher  than  that 

of  uj^.  The  inverse  transform  of  a  function  of  this  type  is 


In  this  equation  all  of  the  roots  X  of  =  0  are  assumed  to  he 
distinct.  This  assumption  is  valid  for  hut  for  and 

vx,  =  0  has  two  zero  roots.  (See  equations  (B2),  (B3),  and  (Bi)-).)  The 
terms  in  the  equations  for  0  and  ijr  resulting  from  the  two  zero  roots 
are 

|g(o)  +  n(o)cT  (b6) 


where 


2 


The  inverse  transforms  of  0^^  from  equations  (B5) 

and  (b6) 


0  ==  +  Age'^^S  +  A3e' 


A^a  +  Ag 


)|f  =  +  B2e'^^^  +  B^e^^^S  +  +  B^c  +  Bg  (b8) 

p  =  +  C2e‘’^2  +  *  C5  (B9) 

The  equations  for  the  rolling  velocity  p  and  the  yawing  velocity  r 
can  he  obtained  from  equations  (B7)  and  (b8)  by  differentiation 

p  =  ^(ApXie'^^l  +  A2X2e'^^2  +  +  A^X^e*^^^  +  A^j  (BIO) 

r  =  +  B2X2e'^^^  +  B^X^e'^^S  +  Bj^Xj^e  ^  +  B^  (Bll) 


where  the  expressions  for  the  coefficients  of  equations  (B7)  to  (Bll) 
are  given  by  equations  (6)  to  (8)  in  the  section  ^Calculation  of  Motions. 
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APPENDIX  C 

SOLUTION  OF  BIQUADRATIC  EQUATION 


Many  methods  are  available,  of  course,  for  solving  for  the  roots 
of  a  biquadratiG  equation.  For  example,  there  are  Horner’s,  Ferrari ’s, 
Bernoulli’s,  Descartes’,  and  Hitchcock’s  methods;  various  methods  of 
solution  by  trial;  and  also  various  graphical  methods  such  as  that 
given  in  reference  1.  Solution  by  trial  in  which  synthetic  division  is 
used,  however,  is  recommended  as  being  the  simplest  method  for  most 
lateral  stability  work.  The  characteristic  equation  for  the  lateral 
motions  of  an  airplane 

+  CX^  +  DX  +  E  =  0 

generally  has  two  real  roots  and  a  pair  of  conjugate  complex  roots.  For 
these  cases  the  two  real  roots  can  be  factored  out  easily  and  the 
remaining  quadratic  solved  for  the  conjugate  complex  roots.  In  the  few 
cases  for  which  all  four  of  the  roots  of  the  characteristic  equation 
are  complex,  Descartes’  method  can  be  used  to  factor  the  biquadratic 
equation  into  two  quadratics.  When  there  are  real  roots,  solution  by 
Descartes’  method  requires  more  time  than  factoring  out  the  real  roots 
singly  and  consequently  is  not  recommended  for  general  use.  These 
methods  of  solution  are  explained  in  the  following  sections. 


Solution  by  Trial  by  Means  of  Synthetic  Division 

Solution  for  real  roots  by  trial  by  means  of  synthetic  division 
consists  of  successive  approximations  of  a  root  and  checking  by  synthetic 
division  until  the  root  is  determined  to  the  desired  degree  of  accuracy. 
This  check  by  synthetic  division  is  based  on  the  fact  that  if  a  is  a 
root  of  a  polynomial  f(x)  then  x  -  a  is  a  factor  of  f(x)  and  con¬ 
sequently  no  remainder  is  left  when  f(x)  is  divided  by  x  -  a. 

The  method  of  solving  the  stability  biquadratic  equation  by  trial 
with  synthetic  division  is  explained  in  three  steps  in  the  following 
sections.  First,  the  rule  for  synthetic  division  and  a  numerical 
example  are  given.  Second,  the  specific  use  of  synthetic  division  for 
factoring  a  biquadratic  is  illustrated  by  a  simplified  example  for 
which  the  roots  are  known.  This  example  shows  how  the  cubic  and  then 
the  quadratic  factors  of  the  biquadratic  are  obtained.  Third,  the  use 
of  synthetic  division  in  extracting  the  roots  of  a  representative 
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characteristic  stability  biquadratic  is  illustrated  with  special  refer¬ 
ence  to  methods  of  making  the  first  approximations  of  the  real  roots. 

Explanation  of  synthetic  division.-  Synthetic  division  is  explained 
in  gimnat.  all  algebra  text  books  but  is  presented  herein  for  the  con¬ 
venience  of  the  reader.  The  rule  for  synthetic  division  may  be  given 
as  follows : 

Assume  that  a  polynomial  in  x  (f(x))  is  to  be  divided  by  x  -  a; 
write  the  coefficients  of  the  polynomial  in  order,  supplying  0  when 
a  coefficient  is  lacking. 

Multiply  a  by  the  first  coefficient,  and  add  (algebraically)  the 
product  to  the  next  coefficient. 

Multiply  this  sum  by  a,  add  to  the  next  coefficient,  and  proceed 
until  all  the  coefficients  are  used.  The  last  sum  is  the  remainder  and 
also  the  value  of  the  polynomial  when  a  is  substituted  for  the 
variable  x. 

For  example,  divide  x^  +  3^^  +  3^  -  x  -  6  by  x  -  3* 


1  +  3+  3  -  1-.  6 


+  3 
1  +  6 


+  l8  +  63  + 

+  21  +  62  +  180 


Use  of  synthetic  division  in  factoring  out  roots.-  The  use  of 
synthetic  division  to  factor  out  two  known  rational  roots  of  a  biquadratic 
equation  is  illustrated  by  the  following  simple  example.  These  two 
rational  roots  represent  the  two  real  roots  of  the  characteristic  sta¬ 
bility  equation  which,  of  course,  are  not  normally  known  but  can  be 
approximated  by  the  method  given  in  the  next  section  of  this  paper. 

One  factor  of  the  biquadratic  is  x  -  1  so  there  is  no  remainder 
when  the  biquadratic  is  divided  by  the  root  1 


1  +  3  +  3-1-6 


+  1  +  4  +  T  +  6 
1  +  h  +  7  +  b  0 


1 


Since  the  remainder  is  0,  x  -  1  is  one  factor  of  the  biquadratic 

equation  and  x^  +  hx^  +  7^  +  ^  another  factor.  Inasmuch  as  a 
cubic  equation  must  have  at  least  one  real  root,  a  second  real  root  of 
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the  biquadratic  equation  can  be  factored  out  of  the  cubic.  For 
example  x  +  2  is  a  factor  so  divide  the  cubic  by  the  root  -2. 


1  +  4  +  7  +  6 


-  2  -  4  -  6 

1  +  2  +  3  0 


-2 


2 

The  factors  of  the  biquadratic  then  are  x  -  1^  x  +  2^  and  x  +  2x  +  3 
The  quadratic  factor  can  be  solved  for  its  roots  by  the  quadratic 
formula.  For  example 


-2  ±  \/4  -  12 

X  =  - 1 -  =  -1  ±  i 


Example  of  application  to  characteristic  equation. -  Reasonably 
accurate  first  approximations  to  the  real  roots  of  the  characteristic 
equation  can  be  obtained  from  simple  formulas.  Successively  closer 
approximations  can  then  be  obtained  by  interpolating  from  the  remainders 
The  following  example  illustrates  the  application  of  this  method  to 
obtaining  the  roots  of  the  stability  biquadratic.  The  biquadratic 

+  10.43x3  +  16.32x2  +  68. 6x  -  9-10  =  0 


is  of  the  form 


AX^  +  BX^  +  CX^  +  DX  +  E  =  0 

Since  the  coefficient  E  is  generally  much  smaller  than  coefficient  D 
in  lateral  stability  work^  one  of  the  real  roots  (usually  the  smaller 
of  the  two)  is  approximately  equal  to  -E/D  or  it  may  be  more  closely 
approximated  by  the  equation 


or  for  the  particular  case 


_  >^9.10 _ 


0.129 
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Approximating  the  root  by  synthetic  division 


1  •+  10.43  +  16.32  +  68.6  -  9.10 

Appr  oximat i on 

+  .13  +  1.36  +  2.3  +  9.10 

.1284  2 

+  .13  +  1.36  +  2.3  +  9.14 

.129  1 

1  +  10.56  +  17.68  +  70.9  + ^04 

1 

1  +  10.56  +  17.68  +  70.9  +  0 

2 

For  this  root,  the  second  approximation  was  determined  by  dividing  the 
coefficient  E  by  the  fourth  sum  from  the  quotient 


-9.10 
■  70.9 

This  procedure  generally  provides  a  good  second  approximation  for  the 
small  real  root. 

The  cubic  equation  obtained  by  setting 

+  10.561,^  +  17.68X  +  70.9 

equal  to  zero  is  of  the  form 

ax3  +  +  cX  +  d  =  0 

In  most  lateral- stability  work^  a  real  root  of  this  equation  will  be 
approximately  equal  to  -b  or  it  may  be  more  closely  approximated  by 
the  equation 


X  = 


+  d 
+  c 


or  for  the  particular  case 
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Approximating  the  root  by  synthetic  division 


1  +  10.56  +  17.68  +  70.9 

Approximation 

-  9.48  -  10.20  -  70.9 

-9.485  6 

-  9.49  -  10.16  -  71.4 

-9.49  5 

-  9.48  -  10.25  -  70.4 

-9.48  4 

-  9.45  -  10.50  -  67,9 

-9.45  3 

-  9.55  -  9.64  -  76.8 

-9.55  2 

-  9.65  -  8.78  -  85.9 

-9.65  1 

1  +  0.91  +  8.90  -  15.0 

1 

1  +  1.01  +  8.04  -  5.9 

2 

1  +  1.11  +  7.18  +  3.0 

3 

1  +  1.08  +  7.43  +  0.5 

4 

1  +  1.07  +  7.52  -  0.5 

5 

1  +  1.075  +  7.48  0 

6 

For  this  large  real  root  there  is  no  simple  method  of  determining  the 
second  approximation  as  there  was  in  the  case  of  the  smaller  real  root. 
The  magnitude  of  the  estimated  root  in  this  case  is  arbitrarily 
increased  or  decreased  slightly  from  the  first  approximation.  From  the 
remainders  determined  from  the  first  two  approximations^  a  fairly  close 
third  approximation  can  then  be  made. 

Factoring  the  quadratic  equation  obtained  by  setting 

+  l.0T5>^  +  7.48 

equal  to  zero  by  use  of  the  quadratic  formula  gives  the  final  tvo  roots 
of  the  biquadratic  equation. 


X  =  - 


1.075  ±\/l.l6  -  29.92 


-0.538  ±  i 


.  [2C 

W-n 
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=  -0.538  ±  2.68i 
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The  roots  of  the  biquadratic  equation  may  be  checked  by  multiplying 
the  four  factors  to  determine  whether  their  product  equals  the  original 
biquadratic 

(X  -  0.1284)(X  +  9.485)(x  +  0.538  +  2.68i)(X  +  0.538  -  2.68i)  =  (x^  +  9.^57X 
l.220)(X^  +  I.07X  +  7.^7)  =  +  l0.i<-3^3  +  i6.32X^  +  68. 6x  -  9-10 

Solution  by  Descartes*  Method 

Descartes’  method  of  solving  a  biquadratic  equation  is  particularly 
useful  for  solving  equations  which  do  not  have  any  real  roots.  This 
method  is  explained  in  most  text  books  on  advanced  algebra  and  theory ^ of 
equations.  In  general^  the  method  consists  of  reducing  the  biquadratic 
equation  to  a  cubic  equation  which  can  be  solved  easily.  One  root  of 
the  cubic  equation  is  used  to  form  two  quadratic  equations  the  roots  of 
which  are  used  to  obtain  the  roots  of  the  biquadratic  equation. 

Method . -  Reduce  the  general  biquadratic  equation 
AX^  +  4-  CX.  +  DX,  +  E  =  0 


to  the  form 


X^  +  bx3  +  cX^  +  dX  +  e  =  0 


by  dividing  by  A. 

Obtain  the  values  of  q,  r,  and  s  from  the  following  equations 


q  =  c 


-  3  b2 


be  1  ,3 

r  =  d  -  -  +  Q  b 


bd  ^  b^c  3 

®  ®  -  T  ^  ir-  ■  ^ 
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and  form  the  equation 


x6  +  I  qx^  q2  _  1  ^^^2  _  ^  ^2  ^  q 

and  solve  this  cubic  equation  in  for  one  of  its  roots  ^  0. 

Solution  by  trial  by  means  of  synthetic  division  is  recommended. 
Determine  the  values  of  Z  and  m  from  the  equation 


I 


4-  2x' 


r 

’  ^ 


0.  P  r 

"  '  2  +  E 

Substitute  the  values  of  Z  and  m  and  the  value  of  x  used  in 
obtaining  Z  and  m  in  the  equations 

2 

y  +  2xy  +  Z  =  0 


2 

y  -  2xy  +  m  =  0 

and  solve  these  quadratic  equations  for  their  roots  y  from  which  the 
roots  of  the  biquadratic  equation  may  be  obtained  from  the  following 
relation: 
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which  may  he  referred  to  as 


where 
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Figure  2.-  System  of  axes  and  angular  relationship  in  flight.  Arrows 
indicate  positive  direction  of  angles,  =  a  -  c. 


R  =  0 ,  oscillatory 
stability  boundary 
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ure  4.-  Lateral-stability  boundaries  calculated  in  table  I.  Ci  vas 

ip 

the  dependent  variable.  was  the  independent  variable.  wa 

actually  varied  by  changing  Cy.  •  Varying  0^0  in  this  manner 

^tail  ^ 

caused  changes  in  the  tail  contribution  to  all  the  other  derivatives. 


Aspect  ratio,  A 


Aspect  ratio,  A 


Figure  5- -  Variation  of  lift-curve  slope  with  aspect  ratio,  taper  ratio 
and  sveepback  for  the  case  of  subsonic  incompressible  flow,  aQ  =  0. 
Values  from  reference  3^- 
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Figure  6.-  Effect  of  horizontal-tall  location  on  the  effective  aspect 
ratio  of  the  vertical  tail  case  of  subsonic  incom¬ 

pressible  flow,  a  =  0°.  Taken  from  reference  35. 
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Figure  9- -  Effect  of  dihedral  angle  on  C7  for  the  case  of  subsonic 

‘p 

incompressible  flov.  Taken  from  reference  58  and  66. 
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Figure  11.-  Charts  and.  formulas  for  estimating  C7  for  the  case  of 

ip 

subsonic  incompressible  flow.  Taken  from  reference  89. 


_ - _  (1  +  2  sin2AA-t2cos_A 

jtA  cos^A  ^  A  +  4  cos  A 
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f^p)cT=0 


A  +  4  cos  A 

CL=0,ao=2jt  /"—V  +  4  cos  A 


where 


Aspect  ratio,  A 

Igure  12.-  Charts  and  formula  for  estimating  Cnj.  for  the  case  of 
subsonic  incompressible  flow.  Taken  from  reference  25* 


Aspect  ratio,  A 


Figure  13.-  Charts  and  formula  for  estimating  C7  for  the  case  of 
subsonic  incompressible  flow.  Taken  from  reference  85. 
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Theory  (references  26^  34,  and  107} 

- Fairing  based  on 

experimental  data 


Figure  l4.-  Examples  of  suggested  fairing  of  theoretical  values  of  lift- 
curve  slope  for  use  in  estimating  values  for  the  vertical  tall  in  the 
transonic  range. 
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